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SYNOPSIS 


HEMMBA KUMAR DEY SARKER 
Ph. D. THESIS 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 


STRESSES AT THE ROOTS OF CRACKS, NOTCHES 
AND CAVITIES, AND MECHANISMS OF BRITTLE 
AND DUCTILE FRACTURE OF METALS 


The present work mainly considers certain meehanlstie 
aspects of brittle and ductile fracture of materials in the 
presence of cracks, notches and cavities. 

Firstly, a method of finite integral transform is 
developed to resolve the boundary conditions for a wedge 
type plate with radiused notch at the vertex. The stress 
intensity factor for a crack, radially emanating from the 
root of a semi-circular edge notched plate in tension, is 
calculated by this method. These SIF results are -ttien verified 
by the transmission photoelasticity experiments on model 
specimens. 

Secondly, the elastic stresses at the root of a 
crack"like-notchciti« found by applying this finite integral 
transfom method. The geomet:^ of this external notdti is 



chosen in such a manner that the notch is deep enoiigh to 
satisfy the Griffith's condition, but at the same time the 
applied stress at the tip may be raised to a high value in 
the order of fractiire strength of the material. The operable 
mechanism of fracture at the root of such a crack-like-notch 
has been studied for elastic as well as elastic-plastic 
materials* 

Thirdly, for an incompressible elastic solid under 
uniaxial tension with small strain, the exterior interaction 
strain fields between two adjacent prolate spheroidal holes 
are determined from the multipole expansion of the displace- 
ment potentials. The shear strain localization in the 
ligament is obtained for a number of prolateness ratios and 
relative angular positions . At a certain orientation, a 
shear buckling type of folding instability is seen to occur 
at a particular point in the ligament. In this critical 
position of voids, not only the shear strain at that folding 
point in the ligament becomes infinitely large, but also 
the stress at a particular point in equatorial plane on the 
void surface increases indefinitely. 

Finally, these results for elastic solution on 
spheroidal holes and crack- like- notches, are utilised to 
understand the mechanisms of brittle and ductile fracture 
of metals, j 

The present analysis shows, when a constant pressure 
opening loading is applied on the flat faces of a cracto-llke 



charpy type notch, a mismatch radial stress on 0°-plane 
as well as a mismatch shear stress on 70°-plane will he 
produced. For a brittle elastic solid, this mismatch 
radial pulling force causes the notch tip to expand. For 
an elastic-plastic material, the mismatch could produce a 
generalised Somigliana dislcxiation on the crack extension 
plane , The d ecrease in misfit energy of the body has 
maxima first at 0.2° and then again at 0.6° notch opening 
angles, Irrespective of the tip blunting. Experimental 
observations on short edge- cracked beams in different 
metals confirm this, showing two distinctly different 
instabilities for elastic- plastic materials. There are 
five distinct types of stress fields at the root of a 
crack- like-notch : a small core region very close to the 
tip followed by a crack-like stress distribution and then, 
three compressive stress regiors at a far field due to the 
bending action. For a generalised plane stress situation, 
the radial stress becomes compressive on 47 i plane 

s 

whereas the tangential component of the stress continues 
to be tensile ; then around 70° plane both the components 
become compressive, leaving a pure shear field over a small 
portion on this plane. 

When this model is extended to an elastic- ide ally- 
plastic Von- Mises material, under plane strain condition, 
a large shear localization is seen to occur on 45° plane, 
and a high plastic constraint is developed below the notch 



root on 0° plane. For this material., in the absence of arQ?' 
work-hardening, the maximur^ plastic constraint factor could 
easily he as high as five or six. The role of elastic 
mismatch stresses iS discussed in relation to the crack 
extension on 0° and 45° planes, depending on whether the 
local instability occurs in a region of high dilatation or 
pure distorsion with hydrostatic tension. Due to the abrupt 
release of misfit energy the first local instability can 
take place in these stagnant flow regions, nucleating micro- 
cracks or micro-voids close to the last vestige of elastic 
cores. The Criterion for discontinuous ductile dimple 
fracture ahead of a crack- like-notch largely depends on the 
stability of growing voids and their coalescence due to 
high shear strain localization in the ligament between 
adoacent spheroidal holes. For a Von-Mises ideally plastic 
material, based on slip-line theory, a simplified analysis 
on the stress strain fields between two spheroidal holes 
reveals that there is a thin layer of material in the 
ligament which does not require any further deformation. 

As there is a sharp shear strain gradient in this region, 
tjjis would be favourable for local plastic necking, causing 
ductile fracture. 

These theoretical results agree favourably with 
the experimental observations. 



CHAPTER 1 


INTRODUCTION 


Strength and toughness are usually regarded a.s 
the two most important mechanical properties to select a 
material for an engineering design. These properties of 
a material are not determined by an average behaviour of 
regular atoms, but mainly by their gross disorder in the 
aggregate. Energy required to create a new surface during 
the fracture of a truely brittle solid is obtained from the 
area under the force-displacement relation between the 
pairs of atoms representing the cohesion of the solid. 

For these materials, the distance between the two cracks 
surfaces at the tip is in the order of an interatomic 
spacing at the onset of a brittle fracture, giving a stress 
level as high as the cohesive strength of the materictl. 

On the contrary, it is observed that the TMaximum attainable 
stress at the root of a machined notch is well below the 
material cohesive strength. Thus, there is a need for a 
careful study of a crack like notch problem from the mechanics 
of elastic continuum. 

Pig, 1,1 shows the geometry for such a crack 
like notch. The flank angle and tip-root radius are small, 
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but, always finite,. This geometry is very similar to a 
charpy type notch near the root region, but it is deep 
enoi^h to satisfy the Griffith's condition ^1, 23 at the 
onset of fracture, 

1.1 Griffith Criterion for Brittle Fracture of Elastic 
Solids 

Griffith first made a systematic study on the 
problem of brittle fracture of solids in the presence of a 
crack. The Griff ith-equation Cl, 23 is based on an 
energy balance criterion. So, Griffith considered the 
First Law of thermodynamics at the onset of crack instabi- 
lity leading to catastrophic fracture. Fig. 1.2 shows the 
geometry considered by Griffith. Accordiig to Griffith 
[Il» 23 j the crack will propagate under a constant applied 
stress if an incremental increase in crack length produces 
no change in the total energy of the system. This means, 
the increased surface energy is compensated by a decx’ease 
in elastic strain energy at the onset of crack instability. 
Thus, a crack will propagate when the decrease in elastic 
strain energy is at least equal to the energy required to 
create the new surface. Fig. 1,2 shows a centrally through- 
cracked plate having crack length 2G which is loaded by a 
remotely applied stress a , For this geometry, the Griffith *s 
condition at the onset of fracture, gives, <^o “ where 
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~ ^ plane stress (1.1) 

2 Ey 

=/ |- fQj. plane strain (1.2) 

ttGCI - v^) 

where Y is the siipface energy per mit area, E and 
are Young’s modulus ani Poisson’s ratio of the material 
respectively. 

Griffith’s equation is a necessary condition 
for fracture, hut not always sufficieit for the growth of 
a crack Csll • The additional condition to be met is that, 
there must be an operable mechanism of fracture at the tip 
of a crack. In our crack- like-notch problem, the geometry 
has been chosen in such a manner that the notch is deep 
enough to satisfy the Griffith’s condition, but at the same 

t 

time due to the tip bluntirg , the applied stress at the tip 
may be raised to the fracture strength of the material. 

The problem of wedge opening of a crack which 
is explained in Fig. 1,3 has been studied by several 
investigators L3 - 83 If a crack is sprung open by a 
fixed wedge Inserted between its faces, the equilibrium is 
stable, because the wedge applies a load on the crack face 
in a direction normal to i^e wedge, and hence, as the wedge 
is moved in, the wedge load poiit does not suffer any 
relative displacement in the direction of the load. The 
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crack grows to a fixed length governed hy the thickness of 
the insert. There would be a misfit energy created by the 
fixed insert of given thickness on the plane of crack 
extension which must be compensated by the generation of 
new surfaces at the crack-tip. Measurement of stable crack 
length by this wedge openirg method provides a means for 
determining the surface energy of brittle materials* 

1.2 Barenblatt*s Hypothesis at the Tip of an Equilibrium 
Crack in a Brittle Solid 

The equilibrium and finiteness of stresses at 
the tip of a crack wefe introduced by Barenblatt 
According to Barenblatt, the tensile stress at the contour 
of an equilibrium crack is finite, and the opposite faces 
of a crack close smoothly at the end where forces of 
attraction (or cohesive forces), are present in small zones 
(Fig, 1,4 (a) ). These cohesive forces pull the crack faces 
together and form a smoothly joined cusp of two crack faces. 
Griffith was aware of this, and mentioned about the "molecular 
attractions" across such a crack very close to the tip. 

The width of this cohesive edge region of a crack is very 
small compared to the size of the whole crack. The form of 
the normal section of the crack surface in this cohesive 
edge region does not depend on the acting loads, and is 
always the same for a given material under given conditions. 
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The forces of cohesion that act at the surface of a crackj 
compensate the applied extensional loads and secure finite- 
ness of stresses with the smooth closing of the crack faces. 
With 8 J 1 increase in extensional loads, forces of cohesion 
grow, thus s-djusting themselves to the increase in tensile 
stresses. The crack starts expanding only upon reaching 
the highest possible forces of cohesion at the edge. The 
successive e^q^ansions of the crack edge under increasing 
loads are schematically shown in Fig. 1,4 (b). It will be 
shown later, that if a crack- like-notch does not expand, 
and the plane section remains plane near the tip, then a 
generalised Somigliana dislocation will be produced on the 
crack extension plane, 

1,3 Elastic Stresses at the Roots of Cracks and Notches 
Fig. 1,5 shows a schematic representation of 
the stresses on zero degree plane for a notch as well as 
a crack for comparison purpose. The infinite stress at 
the tip of a crack has arisen due to the assumption that 
the crack tip is very sharp, and the tip root radius is 
exactly zero. The inverse square root singularity of 
stresses would be valid near* the tip region CSO - 27 3* 

But, actually the tip edge boundary is always traction free. 
Hence, there is a need to study the problem of a crack-like 
edge notch where the crack is infinitely deep \d.th some 
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finite smell flenk °ngle end tip-root rediue blunting . 

Elasticity equations consist of the equilibrium 
equations, compatibility condition and the boundary 
conditions. The very basic fundamental assumption for 
deriving the equilibrium equation is that, the stresses 
and strains must remain finite everywhere. Thus, the 
existence of a strong discontinuity in the formulation of 
the equilibrium equation is ruled out. However, in the 
theory of linear homogeneous isotropic Hookeaii elastic 
solid, a weak discontinuity in the form of a "generalised 
Somigliana dislocation" is well recognised in the frame 
work of infinitesimal strain thsory - 13 2]. Displace- 

ment^^umps like kinks, foldings etc. are well Imown in the 
theory of elasticity. Furthermore, excluding the possibi- 
lity of generating discontinuity in the form of shock waves 
or dislocations, it is an impossibility to have a stationary 
discontinuity in an elastic solid. The second question 
arises whether a weak stationary discontinuity is possible 
on the crack extensional plane in an elastic solid if the 
gradient of stress goes to infinity at a point with conti- 
nuous finite stresses everywhere. Thomas C 12 3 and Hill 
Ci3 3 have shown that such local abrupt strain oumps 
necessarily propagated. 


are 
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1.4 Brittle Practure at the Root of a Crack-Like-Notch 
in an Elastic Solid 

Ohreimov^s cracked body wedge opening 

problem is of great practical importance in fracture 
mechanics [l29 3» Cottrell CsII equated the misfit 
energy created by a fixed insert of thickness b^ during 
wedge opening of a crack, ^^rith surface energy for stable 
crack growth. The length of the stable crack growth can 
be calculated as follows : 

2 

Misfit energy = G b^ / 4Tr(l - v ) 

Sijrface energy = 2 Y C 

Hence, the stable crack size is given by, 

2 

C = Gb^ /Sud - V )Y (1.3) 

where C = equilibrium length of the stable crack extension 

Y == the energy required to create a new surface 

V = Poisson *s ratio of the material 

G = the shear modulus of the material. 

t of stable. Cr^cK ^rooxtP{ 

The essehtial features of equilibrium stable 

crack growth during wedging are incorporated while formula- 

lilct’- 

ting our crack-notch problem. The present work ccosiders 
only a finite radi used- tip notch with small flank angle. 

The author considers this as a more realistic geometrical 
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configuration for a crack. This probl-m is analogous to 
hydraulic fracture of rocics and oil-bearing stratum Cso]]* 

Thus, there must be tvo distinctly different 
conditions for fracture of a crack like notch Csi 3. 

Firstly, the mechanistic equilibrium condition satisfj’'ing 
the finiteness and continuity of stres-es, strains and 
strain- derivatives everywhere, and secondly'-, thermodynafaic 
condition for an energy balance at the crack- tip during the 
onset of fracture. In man^r problems the Griffith condition 
may be sufficient. However, for a crack-like -notch, if 
one satisfies all the conditions of stresses and strains, 
the mechanistic condition causes either a radial expansion 
of the notch tip (or produces a Somigliana dislocation on 
the crack extension plane due to a mismatch radial stress). 

If this expansion is taken as a small crack, emanating from 
the radiused notch tip, like Bowie’s problem Z322, a 
brittle fracture criterion for a crack- like -notch can be 
proposed. Alternatively, one could think of a force of 
attraction on this dislocation due to its image Css]]]* 

1,5 Bowie’s Problem and Stress Intensity Factor for Cra,cks 
Emanating from an Edge Notch 

In an actual service situation, small cracks are 
generally seen ahead of stress concentrations like notch-roots, 
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key- ways etc. Failures are very often seen to occur due to 
the initiation of a small crack emanating from these machined 
notch-groove stress raisers. As soon as a very small crack 
in the order of few grains starts from a notch root surface, 
stress field pattern in the root region alters drastically. 
Bowie [132 3 consid!=red the fracture behaviour of crocks 
emanating radially from the diameter of the circular hole 
in a tension plate as shown in Fig. 1.7. It is very important 
to understand the effect of notch geometry on the stress** 
intensity factor for such cracks. The participation of 
notch talcing an equivalent crack length should be known for 
such geometries [[34 - 37 3* 

Fig, 1.8 shows the geometries and the loading 
patterns of the problems considered. The S.I.F. of such 
cracks has been calciilatsd by using a finite Mellin transform 
method *" 18 3) and then reducing the problem to a 

Weiner- Hopf integral equation, which is finally solved by 
factorization principle [[383 • These results can easily be 
verified by transmission photoelasticity experiments. In a 
brittle material, at the onset of fracture, the notch tip 
starts expanding for a crack like notch, and the fracture load 
can be calculated ly S.I.F. methods used in Bowie *s problem. 
Xhese ealculations are found to be very useful to predict the 
fracture load for a crack in a , notched body , 
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1,6 Elastic-Plastic Stresses Ahaad of Cracks and Notches 
Under A Plane Strain uondition 

In metalsj plasticity- spread at the root ' 
of a crack is well estahliehed 39 - 44 3. For high 
strength raaterialsj or raaterials with small , scale yielding 
property, the plastically yielded zone is normally small, , 
and coTild he neglected. From Irwin C4S 3 and Orowan’s 
observations, it is kno\<jn that in Griffith's equation the 
constant ( = '^a ) is a function of the plastic energy 

dissipation rate and the energy . required to create a new 
surface at the onset of fracture. However, in metals, energy 
required to create a new surface can be ignored and the 
Griffith's constant which is a material property can be taken 
as entirely due to the plastic en<5rgy dissipation rate. _,¥hen 
the plastic zone is very small compared to. the crack size, 
then erther"^ a simple LEFM,, approach, or Irwin-McClintock^s 
plasticity correction [142 3 on crack length, may be suffici- 
ent to predict the fracture load. Some tip blunting, flank 
angle opening and plasticity-sprepd below the root are always 
noticed in a metal, before the crack becomes unstable in 
static fracture. So, in an elastic- plaetic material it 
would be more appropriate to treat the geometry of a crack- 
tip as a blmted notch, rather than always, a sharp tip, 

A brittle fracture ■ ’of mild steel plate can take place under 
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a plane strain situation, and the local tensile stress of a 
metal ahead of a crack tip <i3Xi be several times the yield 
strength of the material Csi? “■ During fracture 

tests of a short cracked beam , it is found that, it would 
be better to treat the geometry as a crack“like-notch L5 i 3. 
Recent experiments further confirm the notch-like behaviour 
of a short crack during the fracture of elastic-plastic 
materials - SS]]] • Fig* 1*9 shows a typical load-COD 

diagram and crack geometry observed at the tip together with 
the plasticity spread during the fracture tests of a short 
cracked steel plate. In metals, during the plane strain 
fracture of a cracked beam, the plastic constraint and shear 
localization inside the plastic zone normally play the most 
important roles 54 , 55 []]. One must also distinguish 

between the yield stress for general yielding and that for 
local yielding in a region of concentrated stress. 

Experiments C56, 583 on mild steel notched bars at low 
temperature have shown a brittle cleavage fracture which is 
preceded by deformation twinning, and a deformation twinning 
is further preceded by a slip in the grains ahead of the 
notch . 

For an ideally plastic material, under plane 
strain deformation, the maximrcn stress elevation below a notch 
root depends on 359 - 64 3 • yield strength of the 

material, notch flank angle, and the tip root radius. The 
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FEM approach adopted hy the previous investigators to find 
the elastic- plastic stress fields ahead of a notch (Pigs* 
1*10 and 1,11) is far from complete. The major difficulty 
in applying FEM analysis for such a crack like notch 
geometry arises due to the matching of the elastic far 
field with the "bifurcations of plane strain plasticity 
equations. A pure distortional stress in an elliptic field 
could arise, demanding a bifurcation of hyperbolic differ- 
ential equation of plane strain plasticity inside the 
yielded zone £ 502 . It is kno\*jn that, during elastic- 
plastic torsion of a prismatical bar, there could be an 
elastic core inside the plastically yielded material 
|“50, 65 - 673 • Hence, in strict formulation, the problem 
is extremely complex, and possibly, the characteristics 
method would only give an accurate elastic-plastic stress 
field, "Ghder plane strain deformation, a pure shear field 
with hydrostatic tension on 46*^ plane could exist, for a 
large scale yielding ideally plastic Mlses material 

Css, 123 3* 

1,7 Mfechanism of Elastic-Plastic Fracture at t he Root of 
a Crack-Like**Notch 

The purpose of measuring the fracture 
toughness of a material is to determine a single parameter 
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characterisl.g its resistance to fast fracture . The 

crack toughness is regarded as the critical value of the 
strain energy release rate for that material at the onset 
of fracture. Thus, the toughness is also a measure of the 
resistance to crack extension. The load at which a crack 
extends needs to be determined with reasonable accurac5% 

Irwin- McClintock’s plasticity correction 
method, ASTM recommended secant shift procedure on LSPM and 
Rice’s [1693 path independent J-integral in terms of pseudo- 
potential energy have been partially successful for small 
scale yielding materials only. 

Thus, it is necessary to understand the 
nature of elastic as well as elastic- pla.s tic stresses ahead 
of a crack- like- notch. So, in an elastic-plastic solid, the 
plasticity spread is neither all that small, nor the effect 
of tip blunting can be ignored. When the fracture stress 
approaches the yield strength of tiie material, the effects 
of tip blunting , plasticity spread and notch toughness 
rather determine the fracture behaviour. Brittle fracture 
in steel has, however, bean observed at stresses well below 
the yield strength of the material. 

^ Experiments conducted on charpy 'type edge 
notched "beams give an idea about the overall deformation 
produced by localized slip. Using Griffiths-CXren £593 
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results on charpy notch, Knott Smith, Cook and 

Ran C Z explained the nature of cleavage fracture 

at the root of a crack. Bilhy, Cottrell and Swinden 
calculated the plastic strain aJiead of a sharp notch from 
dislocation theories. Experimental observation ahead of 
yielded notch shows, micro-cracks, micro-voids, flow- 
localization, glide band decohesion and other types of 
localized deformations in the regions of high dila.tation 
and pure distortion. Under plane strain situation, in 
some materials micro- cracks are produced near the elastic- 
plastic interface on 0° plane, wheraas in other materials 
micro-voids and their coalescence are seen to occur first 
on 45° plane inside the plastic zone [[[72 - 73 3* Metal 
intrusions and extrusions are also visible on the free 
surface at the onset of crack growth. 

It will be shown later, that the nature of 
mismatch stresses observed at the root of a crack-like- 
notch is in the form of strain derivative jump. This type 
of weak discontinuity can travel freely like waves and has 
a tendency to produce a displacement jui-^ip and velocity 
discontinuity at the elastic-plastic interface 
This mismatch cen cause asymmetry of the loading as well. 
As t^ere are t\-ro stagnant flow separ=^tion regions in a 
ful|-y developed plastic flow pattern ahead of a crack- like 
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notch, velocity discontinuities in plastic flow ^^3 

can occur in these regions where there are sharp changes 
in the boundary from elastic to plastic . The 

discontinuity in chear strain-derivatives c^n also occur 
in these regions 1112, 13, 743* 

The discontinuous plastic yielding of 
steel produces a sudden large avalanche of dislocations, 
so that the local stress can build up very rapidly which 
might create cleavage micro-cracks ahead of the notch 
3 76 - 783* There must be a driving force to nucleate a 
micro- crack in the high dilatational region aliead of the 
notch, and then this micro-crack needs some applied stress 
for its groirth. Thus, under plane strain condition, in an 
elastic -plastic materia3-,the crack extension can take place 
either on 0°plane due to the maximum normal stress in a 
hyperbolic-parabolic field, or on 45® plane due to the 
shear stress component elevation in a parabolic-elliptic 
field. 

Fig. 1.12 sho\-js the nature of the formation 
of micro-cracks and micro- voids inside the plastic zone 
ahead of a crack- like notch. 
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1,8 Interaction Strain Field Between Spheroidal Cavities 
And Mechanisms of Ductile Rupture Due to Shear 
Localization 

In ductile metals, fracture normally occurs 
by the growth and coalescence of holes. During tensile 
test of a round mild steel bar, small three dimensional 
holes are cotnraonly seen to occur in tlie necked region before 
final fracture. Orowan first pointed out that at 

the centre of a round specimen, there would be a large quadri- 
lateral shape prolate spheroidal void in tiie middle portion 
of the necked region. This can be explained as follows s 
from tensor algebra, a uniaxial tension can be resolved 
into a hydrostatic tension plus two 45^ shaars Thus, 

the hydrostatic tension at the centre of a necked tensile 
specimen helps to create void eround second phase particle, 
and the rotation ally syraraetric 45° pure shear stress deforms 
the void into a spheroid. At the onset of void instability, 
the prolate-spheroid becomes a quadrilateral shape CsiU . 

In recent years, considerable attention has been focused on 
the stability of .a single growing void and the ligament- 
fracture-instability of an array of voids in plastically 
yielded material The growth and stability of a 

three dimensional spheroidal hole will be different from the 
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plane growth of cylindrical holes. A cylindrical hole in 
a plana strain situation, can grow very rapidly beyond a 
certain nominal strain giving a quadrilateral shape Cain. 
A single prolate spheroidal hole in an incompressible 
elastic solid can also become completely unstable violating 
the energy balance criterion, especially when the prolate- 
ness ratio reaches a particular value []] 823 • 

McGlintock [Issj and Berg 11^4 3 studied 
the plastic flow and shear localization between two cylin- 
drical holes considering the traction-displacement boundary 
conditions for plastic fracture 3353* Ihe plastic flow 
localization at a point would occur when a non-deforming 
surface through that point, and some field of relative 
motion across a liiin layer parallel to that surface, does 
not require any further deformation in the material on 
either side of the layer (Fig. 1.13). 

As elastic stresses in an incompressible 
solid give a prior history to the elastic-plastic fields 
at the incipience of yielding, it might be worthwhile to 
study the exterior interaction strain field between two 
adjacent prolate spheroidal holes in tension. Pig. 1,14 
describes the geometry and tension-loadir^ for a pair of 
prolate holes in space inside an incompressible solid* 



18 


1,9 Objective of the Present Work 

The purpose of undertaking this work is, 
firstly* to study the brittle as well as elastic-plastic 
fracture behaviour for a crack-like -notch geometry, and 
secondly, to understand the mechanism of ductile fracture 
due to the shear localization between two adjacent holes. 

So, in the present thesis an attempt has been made to 
investigate the mechanisms of brittle, quasi-*brittle and 
ductile fractures purely from continuum mechanics. However, 
there are limitations in applying continuum mechanics on 
such problems £882. Thus, this work is merely an attempt 
to explain the brittle fracture postulates of Griffith* 
Barenblatt, the elastic-plastic fracture concepts of Irwtnr 
Orowan-Bilby-Cottrell, and the ductile fracture theories 
proposed by Berg-McGlintock. 

Firstly, the advantage of using finite 
Mellin transform is demonstrated, mainly to resolve the 
boundary conditions for a notched plate more easily. Then 
a general formulation of elastic stress field at the root 
of an external deep crack like notch is presented. The 
effect of notch angle on mismatch stresses and misfit 
energy is emphasized. 

Secondly, the stress tn-tensity factors for 
a crack emanating from a eemi**circular edge notch are given. 
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Analytical S.I.F. result obtained are ’'.’’erified by 
transmission photoelasticiiy experiments. This method 
is uhen used to find the load carrying capacity of a 
crack-like-notched plate. 

Thirdly, for a Von-Mises ideally plastic 
material, the elastic-plastic stress field at the root 
of a crack-like-notch is presented, showing the regions 
of high dilatation and pure distortion. The roles of 
mismatch stresses are explained. The nature of velocity 
discontinuity and local plastic instability in relation 
to the crack extension is also discussed. 

Finally, the theory of ductile fracture 
Is studied considering the ligament instability of two 
adjacent voids. A highly simplified model is proposed 
for the shear localization between two prolate spheroidal 
holes. The entire analysis is based on the exterior 
interaction strain field in an incompressible elastic 
medium. This model is then extended to explain the 
velocity discontinuity and local plastic instability. 

Void surface critical normal stra3ja is also examined for 
a number of prolateness ratios. 

The above theoretical results On crack-* like-* 
notch and spheTOldal caVlty interactions are compared with 
experimental observations of others. ' 



CHAPTER II 


LITERATURE REVIEW 

2.1 Theories of Brittle Fracture 

A hrittle fracture is normally characterised hy a 
catastrophic failure associated with fast propagation of 
cracks* The plastic deformation at the tips of cracks or 
notches is either negligibly small or practically zero. 
Normally I in high strength materials such types of fracture 
are commonly seen, Griffith Cl, 2^1 first considered the 
energy balance at the onset of fracture. Griffith was 
fully aware that this is a necessary condition, but not 
sufficient, because at crack- tip the equilibrium condition 
should also be satisfied. To satisfy the equilibrium 
condition, Griffith made an attempt to improve the descrip- 
tion of the crack-tip, and according to his own estimate 
the magnitude of the radius of curvature with a finite 
radiused tip, at the end of the crack is in the order of 
intermolecular distance. 

Mott gave an excellent discussion on the 

real atomic structure of a crack- tip. Barenblatt 1^93 

first investigated the equilibrium condition and on 
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crack-tip shape systematically from both continuum and 
atomistic points of view. Barenblatt pointed out, 

if the equilibriiim equations of elasticity were to be 
satisfied at the crack- tip maintaining the finiteness of 
stresses, strains and displacements, then the tip of the 
crack must expand (Fig. 1.4(b)). Fig. 1.4(a) shows the 
shape of Barenblatt’ s equilibrium- crack which is different 
from Griffith’s elliptic crack of finite tip root radius. 

The opposite faces of the crack must close smoothly and the 
tensile stress at the tip of a crack is always finite. The 
rounded tip must expand and form a cusp at the tip to 
satisfy the equilibrium conditions together with the 
finiteness of stresses, strains and displacements. 

Obreimov j Mott tl493, bailey 3 and Benbow and 

Roesler considered the stable crack growth by 

inserting a finite thickness wedge which provided a means 
of determining the surface energy of brittle materials 
such as mica. Cottrell Zsl emphasized that in this 
method, the misfit energy created by the fixed Insert was 
expended in generating a new surface for stable crack growth. 
Barenblatt an<3. Frenkel suggested that the 

mismatch energy on the crack extension plane would be 
responsible, either for the expansion of the crack-tip, 
or for generation of the elastic mismatch displacement. 
Recently, Cherepanov CssQ, Cook, Rau and Smith [90 
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studied the elastic raismatch ahead of a crack- tip as well 
as oi- the interface of tvo-phase materials. Cottrell CsO, 
Bilhy, Cottrell and Swinden Bilby and Eshelby 

QllJ represented the crack-tip deformation by an array 
of dislocations. Again, it is known from Eshelby *s CloH 
analysis on mismatch in an elastic medium that a generalised 
Somigliana dislocation has a displacement discontinuity 
with continuous stresses and strains. This means that the 
stresses and strains are finite and continuous everywhere, 
the strain- derivative may be discontinuous at a point • 

Hill [I132], Bilby [1^1 Us Thomas and Cherepanov 

[^89 3 gave considerable emphasis on such folding type of 
‘ displacement jump in the theory of fracture instability of 
solids. 

Our analysis on the elastic stress field of a crack- 
like axternal deep notch shows that such a generalised 
Somigliana dislocation can be generated on the plane of 
crack extension. This means, there will be a misfit radial 
stress developed on the crack extension plane pulling the 
crack>*tip outward for its expansion. This misfit is 
obtained entirely due to the geometry of a cracls-like-notch 
upon loading. As the notch flank angle is increased, this 
mismatch on the cracb-extension plane decreases rapidly. 

It will be shown later that the creation of this type of 
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mismatch stress due to loading has to he the operable 
mechanism, in addition to satisfying the Griffith's 
energy balance criterion at the onset of fracture. 

2.2 Stress Intensity Factor for Cracks Emanating from 

Machined Notches 

Machined notches, grooves and undsr~cuts are practi~ 
cally unavoidable in designing machine members. Such 
type of abrupt change in the geometrical configuration 
always produces stress concentration. The majority of 
service failures due to fatigue-cracks normally originates 
from some form of notches or stress concentration regions. 
The extent of the notch stress field influencing the crack 
growth behaviour is of great practical importance. As 
the basic parameter in linear elastic fracture mechanics 
is the stress intensity factor for a crack, one must 
evaluate the influence of notch stress field on S.I.F. 
for such cracked bodies. The initiation of this type of 
fatigue-crack can be considered in two stages. Firstly, 
a definite number of cycles is required to fatigue-harden 
the notch root, and secondly, a small visible crack is 
initiated from the machined boundary of the notch. Finally, 
this minutely small crack propagates under the application 
of the notch stress field, until it becomes so large that 


■I 
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the effect of notch geometry on S.I.P. can he neglected 
completely. 

Winnie and Wimdt E 92]] emphasa'’ed the practical 
importance of this problem in the case of failure for 
turbine hubs . The failure analysis of Hinlcly~ Point E 93 E 
turbine revealed a short fatigue crack propagated from the 
root of a key-way notch of a turbine hub. Bowie CssH 
analytically obtained the stress intensity factor of cracks 
emanating from a hole along the transverse diameter in a 
tension plate . Newman C342 examined the nature of S.I.F. 
variation for cracks emanating from various holes and cut- 
outs. Isida and Nisitani C94]]> and Yokobori, Kamel 
and Konosu E95 3 considered the geometry of cracks 
emanating from elliptic holes and determined the stress 
intensity factors in Mode- I and Mode- II I respectively. 

In a recent compendium, Rceke and Cartwright E 96 ]] 
reported many such results for the S.I.F. of cracks 
emanating from various types of cut-outs and holes. Smith, 
Jerram and Miller determined the S.I.F, of cracks 

radially emanating from the root of edge-grooved machined 
notch in a pin-loaded tension specimen. The participation 
of notch geometry with fatigue crack length had been 
inve stig ate d by Miller C3S3 considering an equivalent 
crack length. According to Smith and Miller [1973 at the 
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instant of initiation of a crack from a notch root, the 
participation of notch depth to crack length could he 
zero. However, when the crack grows some distance into 
the hulk material, the equivalent crack length will involve 
both notch geometry and crack size. Broek studied 

the propagation of fatigue crack emanating from a semi- 
circular edge notch in finite plates. His results agree 
favourably with our calculations. Bowie and Neal ClOOU 
and Kutter also considered the stresses at the 

roots of radiused tip • edge slots and cracks emanating 
from various grooves and cut-outs. In all these results, 
as expected, at a certain value of C/R ratio the effect 
of notch geometry on S.I.F. vanishes, the S.I.F. (Kj) is 
then simply calculated as follows: 

For centrally slit cracks Kj = \firX/2 , which is 
obtained from Griffith's result j^l, 2[[]j whereas, for 
an edge crack Kj = 1.1215 \/1 tT , as derived by 

Wigglesworth £"24 3} Koiter Sneddon and 

Sih and Liebowitz E26£l. 

The influence of notch geometry on the size of 
plastic zone at the tip of a crack symmetrically emanating 
from a circular hole, was considered by Sumpter and Turner 
E 1023 . Mubeen and Bandyopadhyay Css: experimentally 
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observed the effects of notch root radius and crack size 
on the shape of plastic zone for a crack emanating from 
a charpy type edge notch. 

Therefore, to predict the fracture behaviour of 
a beam containing a short crack emanating from an edge 
notch, it is highly importejit to determine the S.I.F, 
accurately. Fig. 1.8 shows the details of the geometries 
considered. In the first problem there is a semi-circular 
edge-notch in a semi- infinite plate with a crack radially 
emanating from the notch root but perpendicular to the 
straight edge, whereas the geometry of the second problem 
being a very short crack emanating from the radiused tip 
root of an external deep crack- like- notch. 

2.3 Photoelastic Determination of S.I.F, for a Crack 

One of the easiest ways of direct observation of 
the stress field near a crack-tip is by transmission photo- 
elasticity experiments. Photoelastic technique is a 
powerful experimental tool to study the crack- tip stress- 
strain fields and for the determination of stress intensity 
factors for cracks. Wells and Post [-103 3 and Pratt 
and stock [“1043 studied the dynamic stress distribution 
surrounding a running crack. Wells and Post found that 
the stress fields around stationary and dynamic cracks 
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■were identical. Kblayashi, Harris and Engstrom ClOSH 
determined transient displacements and strain distributions 
in the vicinity of running cracb-tips in centrally notched 
magnesium plates by moire- fringe technique, Dixon and 
Visser used photoelastic coating technique to 

determine the plastic zone size around crack-tip in a 
'centrally cracked tension specimen. Their results agreed 
favourably ■with those of Dugdale under similar 

situations, Gerber ich cio7n measured the thickness 
variation ■within the plastic zone in a centrally cracked 
thin tensile plate specimen. He was, thus, able to deter- 
mine the third directional principal strain from the 
condition of incompressibility and hence, he was able to 
separate all the three principal strains in the plastic 
enclave surro^unding the crack- tip, Kobayashi 
Theocaris IIl09[3j Theocaris and Gdouts 
Vitvitsld-i et al CmH calculated the S.I.F. and also 
studied the elastic- plastic stress distributions surrounding 
a crack-tip* 

Fig. 2,1 shows the coordinate system adopted by 
Williams C 22 I] for which the stresses at the root of a 
crack-tip are given by, 

" vfe* I I 


« • • ' • 
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Irwin foimd an excellent agreement 

"between the above theoretical results and those of the 
photoelastic experiments conducted by Wells and Post [I 103 3 • 
Irwin modified the above equations by adding one more 
constant term to the R.H.S. of the expression for a . 

Raju C 1133 used Irwin *s modified equations to calculate 
the stable crack growth in an elastic- plastic solid. 


From the stress-optic law in a photoelastic 
material isochromatic fringe order H is 

related to the difference of principal stresses ( - ^22^ 

~ ^ "^max » ’'^max maximum inplane shear 

stress, and therefore, 

Vo 

r- 2 ' 2 ' Iffo 

^ ''^.max " n c '^yy ~ '^XX^ ^ II " (2,2') 

where f^^ is the well-known stress-optic constant and 
t is the thickness of the plate. 
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Hence, from equation., (2.1) and (2.2), 




SinG 


(2.3) 


When a plastic coating sheet is cemented on the surface 
of a cracked specimen and loaded, the surface strain 
would he the same as in the plastic sheet. Hence, for 
photoelastic coating technique, 


K, = 


N f E 

2 (l+v) SinG 


(2.4) 


where the fringe value of the coating sheet f is related 
to fg , the fringe constant of the plastic in terms of 
strains hy the expression, f - f g /t^ , v is the Poisson *s 
ratio and E is the Young ^s modulus of the metal prototype. 
Eqns. (2.3) and (2,4) suggest a method of determining 
the stress intensity factor (K^) for a crack hy making 
observations of isochromatic fringe orders at different 
values of r and G in the close vicinity of the crack- tip » 

The accuracy of the photoelastic technique depends, 
on the sharpness of the crack-tip and how closely one could 
measure the fringe orders at the vicinity of the crack-tip. 
Knott [l63[3 gave the percentage error involved in the 
S.I.P. calculations against the distance from the crack- 
tip. Gerherich 107 3 used the modified form of Eqn .(2.1) 
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iDy substituting Kj -o^xfira. where is the applied 
gross stress; and he obtained the following results; 

N ®11 " ^^22 x/''a/2T . Sin 0 (2.5) 

1 ' - f E 

Photoelastic techniques have been applied in the 
present work to study the influences of notch geometry and 
crack length on the stress intensity factor of a crack 
radially emanating from the notch-root. These experimental 
results are then used to verify the theoretical solutions. 

2,4 Elastic- Plastic Fracture at the Roots of Cracks 
and Notches 

In metals, plasticity spread ahead of cracks and 
notches is well-known. For a small scale yielding material, 
this yielded zone is normally small compared to a crack 
length and the load-GOD diagram is almost linear under 
plane strain situation. However, in an actual metal, the 
plastic zone size ahead of the tip of a crack need not 
always be too small. Simplified LEPM approach may not 
predict an exact elastic“plastic stress-strain field in 
such a situation. On the other hand, the results of 45^ 
charpy notch cannot be used either. So, it is as much 
important to know the elastic-plastic stress field as the 
cleavage and fibrous micro-mechanisms of fracture at the 
root of a notch [”56 - 593. Ihis information may provide 
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some lasis for evaluating the fracture toughness of a 
metal [l54 3 • 

The geometry of the crack- tip is very crucial 
and the details of the notch tip shape is somewhat contro- 
versial McClintock CUS^ assumed a V- notch 

shape at the tip to investigate the shear localization 
inside the plastic zone on 45° plane. Fig* 2.2 shows 
Berg-McClintock model - 85]]] on micro- voids formation 

on 0° and 4 : 5 ^ planes. Fracture can take place entirely 
due to the nucleationj growth and coalescence of such 
voids in a highly ductile material. To understand the 
nature of elastic- plastic stress field, overall deformation 
produced hy the localized slip, and to distinguish hetween 
the yield stress for general yielding and that for local 
yielding, it might he more appropriate to assume the 
crack- tip as semi-circular. 

FEM calculations of Griffiths and Owen Cs 93 
showed that the maximum stress elevation would occur very 
close to the elastic- plastic boundary, hut inside the 
plastic zone. Fig. 1,11 shows these FEM results for a 
45° flank angle charpy notch. Rice and Johnson [362 , 

Rice and Rosengren Rice and Hutchinson 

and Rice and Mcmee king [l6lll obtained elastic-plastic 
stress fields ahead of a blunt crack based on a similar FEM 
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calculations. Fig. 1.10 shows these results at the tip 
of a Hunt crack. AcGord:..ig to these authors, the strain 
hardening has a pronounced effect on the plastic constraint 
factor. Kbning [^118 [3 and Anders son reported, 

the flank- angle opening of the notch must also he incorpo- 
rated in the geometry along with the tip blunting. In the 
present work, the author makes an attempt to solve the 
crack- like-notch problem (Fig. l.l). Green and Hundy 
Ewing and Griffiths C12 q3j Bichards and Ewing [3 121 3 
and Rice and Johnson 1 ^ 62 ^ observed that Hill*s [ISOQ 
logarithmic spiral slip-lines near the root of a radiused 
notch would become straight line giving an uniform stress 
field at some distance from the blunt notch-tip (Fig. 1.10). 
Thus, the characteristics of the partial differential 
equations governing the plane strain ideal plasticity, 
would alter from hyperbolic to parabolic. Furthermore, 
McClintock^s slip- line field analysis at the root of a 
crack shows the tendency for a pure shear field at 45° 
plane imder a plane strain situation* 

Our analysis on a crack- like-notch shows that 
a pure shear field would be produced at an angle with 
respect to the plane of symmetry. When this model is 
extended to an elastic-plastic material mder a plane 
strain situation, a pure shear field would lie on 45° plane 
at some distance inside the plastically yielded zone. 
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calculations. Fig. 1.10 shows these results at the tip 
of a Hunt crack. According to these authors, the strain 
hardening has a pronounced effect on the plastic constraint 
factor. Koning [^lis]] and Andersson reported, 

the flank- angle opening of the notch must also be incorpo- 
rated in the geometry along with the tip blunting. In the 
present work, the author makes an attempt to solve the 
crack- like- notch problem (Pig, 1.1). Green and Hundy 
Ewing and Griffiths Bichards and Ewing ClSl]]] 

and Rice and Johnson observed that Hill*s [1^011 

logarithmic spiral slip-lines near the root of a radiused 
notch would become straight line giving an uniform stress 
field at some distance from the blunt notch-tip (Fig. 1.10), 
Thus, the characteristics of the partial differential 
equations governing the plane strain ideal plasticity, 
would alter from hyperbolic to parabolic. Purthermora, 
McClintock^s slip- line field analysis at the root of a 
crack shows the tendency for a pure shear field at 45° 
plane under a plane strain situation* 

Our analysis on a crack- like-notch shows that 
a pure shear field would be produced at an angle with 
respect to the plane of symmetry. When this model is 
extended to an elastic- plastic material under a plane 
strain situation, a pure shear field would lie on 45° plane 
at s OIK distance inside the plastically yielde'd zone. 
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Sokolovsl?y Pi’ag -.r and Hodge 

Nadai C H Hill E [3 established that there would 

be an elastic core inside the plastic region, as a last 
vestige, during the elastic- plastic torsion of a prismati- 
cal bar. Hence, for a plane strain situation, there woiLl-d 
be a tendency to form a hinge-point under bending action, 
inside the plastic zone on 45° plane. Thus the characte- 
ristic root [12231 of the plane-strain (ideal) plasticity 
equations are changing inside the plastic region from 
hyperbolic to parabolic on 0° plane and hyperbolic to 
parabolic and then elliptic on 45° plane, to match the 
elastic stresses on the yielded boundary. Therefore, a 
simple FEM analysis would be inadequate to solve such 
crack- like- notch problems. It would be more appropriate to 
formialate the problem analytically, taking into accomt 
the bifurcation, velocity discontinuity and local plastic 
instability. Thomason presented characteristics 

solutions for hyperbolic-parabolic stress fields at the 
root of a crack under plane stress situation. 

Pig. 2.3 shows Singh *s observations on 

crack- tip deformations and plasticity spread at the root of 
a short edge-cracked beam. Becent experimental observation 
shows, in a ductile material under plane strain situation, a 
crack can extend either on 0°- plane or on 45°- plane (Begley 
and Landes Cvs;: jRice and Sorensen E 733 ). Our analysis 
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proves, under a plane etrain situation, not -only the 

plastic constraint factor could be in the order of five 

for a nonr* strain hardening Von-Mises yielding material, 

but also only two possible manners in which crack 

extension can take place. If the material has a small 

scale yielding property and is more susceptible to form 

micro~craGks and micro-voids at high dilatational region, 

then the crack would naturally extend along 0°-plane, On 

the contrary, if the material is ductile and more prone 

to nuclelite micro-voids under the action of hydrostatic 

tension with a pure distortional stress, then the crack 

can extend along 45°-plane under plane strain loading. 

In a ductile metal, if the crack extension does not take 
o 

place along 45 plane, then a rapid increase of plastic 
zone size on O^plane on the surface would be visible first, 
before showing any noticeable slow stable crack growth 

Kfouri and Miller [[124 3 , Cook, Rau and 
Smith C 903 and Knott [[63 3 proposed a local Griffith- 
like critical separation energy concept as a criterion for 
continuing fracture for a stable crack growth. In a 
ductile material, it may be inappropriate to employ such 
energy fracture criteria, especially, when the plastic 
strain dominates fracture mechanisms near tip region, such 
as nucleation, growth and coalescence of micro-voids or 
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large shear strain localizations. Bilhy 

McOlintock pointed out, the fracture in an elastic 

-plastic material would depend not only on the high stresses, 

hut also on the strain- localizations. Crrsen and Knott O-^CI j 

Rice and Sorensen ^md, Rice and Tracey [Il263 

emphasised that the critical crack opening displacement, 

6^. would control the fracture at the onset of stable 
tip 

crack growth. Cherepanov also presented a similar 

analysis for the crack-tip in the case of an elastically 
incompressible material. Furthermore, at the onset of 
ductile fracture, the tip opening displacement attains a 
size which could be compai'able to the localized slip and 
fract\are process zone (Hahn, Rosenfield and ICanninen C.273). 
In addition to critical COD method and well-known 

J- integral criteria, Paris, Tada, Zahoor and Ernst [l2^ 
proposed a dimensionless crack growth resistance parameter 
T = E(dJ/dl)/ 0'^ in fully plastic specimens. Paris et al 
suggested from experiments on long edge-cracked bend 
specimens that a T-stability criterion would be consistent 
with the observations of stable and also unstable crack 
growths. For stable crack growth situation in a ductile 
metal, R-curve method ha.s also been found to be very useful 
[ 11293 . 

In the present work, it will be shown that the 
mismatch radial stress on 0° plane and mismatch shear stress 
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on 46° plane must play an important role for nucleation 
of micro-cracks and micro-voids in high dilatational and 
pure distortional regions. In a ductile metal, if no 
such micro-defects are produced, then a large shear flow 
would occur first on the 0°- plane at the end of the 
plastic zone before any slow stable crack growth. Experi- 
ments on elastic-plastic materials show, during a crack 
extension, that the crack always propagates from the 
free notch surface to these high stra-ined regions. Hence, 
the fracture of metals is a problem of surface instabili- 
ty due to the metal intrusions and extrusions 

which are normally visible C^^l^ the notch-root 
free- surface before the onset of stable crack growth. 

2,5 Ductile Fracture Due to Three-Dimensional Hole 
Growth Inside a Plastically Yielded Material 

The general behavioirr of ductile fracture 
usually involves the formation of voids (around inclusions 
or second- phase particles) and the subsequent 

growth of these voids to coalesce in an unstable manner 
, Fig. 2.4 shows the discontinuous modes of 
deformation ahead of a notch root produced by growth and 
coalescence of holes... McClintock 83, 1343 j 

13843, Gottrell 3 13 2 3 and Brown and Embury Lael 
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studied the growth of cylindrical holes and the strain 
required to coalesce these holes in plastically yielded 
materials. These models on ductile fracture are based 
on two dimensional deformation of cylindrical holes. On 
the contrary, during ductile fracture, three dimensional 
spheroidal cavities are commonly seen inside a metal, 
particularly, in a plane strain situation. The criterion 
of ductile fracture largely depends on the stability of a 
growing void as well as the shear localization and ligament 
instability between spheroidal cavities. At present, no 
such three dimensional analysis is available on plastic 
unstable hole growth mechanisms, which give a local 
instability of spheroidal surfaces. Bice and Tracey [l3^, 
Brown C136 3 and Nix and Goods C 1373 considered the 
growth of spherical holes in a plastic material. During 
a plane strain deformation state the holes do not remain 
spherical and would always have a tendency to grow in 
spheroidal shape [1134 3, particularly under uniaxial 
tension. 

To understand this three dimensional spheroidal 
hole growth problem, Singh and Bandyopadhyay C8s3 studied 
the shape dependent instability of a single prolate 
spheroidal hole in a linear elastic incompressible solid 
in pure ^ear and in uniax^ tension loading. In a simple 
shear field, the shear buckling strength is given by the 
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shear modhlus times a shape factor "tiie 

presence of a source of di sturhance , a single prolate 
spheroidal hole surface can become unstable under uniaxial 
tensile field. The sourPe of disturbance could arise in 
many -ways e.g. a dislocation or an inclusion nearby, or 
a slight imperfection of the shape of the spheroid. As 
soon as the applied tensile load reaches a critical value 
the prolate hole will try to rotate at this load to 
increase its rotational shearing resistance, but it cannot, 
due to the alignment of its major axis in the direction 
of the applied tension. Therefore, the prolate hole 
becomes imstable in Biot's Cl 303 internal shear buckling 
mode of deformation. In this situation the principal aijcis 
of strain can only rotate [[[iSS]] , When the applied 
tension exceeds Biot's slide modulus of the surrounding 
materials in the presence of the void, the Saint- Venant 
type of torque vanishes and the local torsional stiffness 
is entirely due to the applied initial stress, Eshelby 
£[1393 McClintock pointed out that such three 

dimensional analysis in an incompressible rubber- like 
elastic solid may help to understand the trend of void 
growth mechanism. 

In a plastically yielded material, the three 
dimensional hole growth mechanism is fairly complex. 

Pending such an exact solution, in the present work the 
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interaction strain field "betwesn two identical prolate 
holes is considered inside an incompressible elastic 
solid imder uniaxial tension field. 

Ductile fracture of metals at high tempera- 
ture involves other mechanismsj mainly due to the diffu- 
sion and gas formation inside the cavities. Incorporating 
these complexities along with the three dimensional surface 
deformation together with interaction stress field inside 
a plastically yielded material, would be a formidable- 
task, One has to rely on, or compare with, the experi- 
mental (and simplified theoretical) results obtained by 
McClintock C1343 , Brown and Embury C H ? 3^^ 

Ashby Cmi 3, Tait and Taplin and El-Sudani and 

Knott 3 • Orowan’s observations at the 

centre of a necked tensile specimen as well as Perra and 
Finnie * s C8i3 exparimental investigation on plane 
strain hole growth, explain the shear' dominant instability 
of voids, Bilby, Eshelby, Kundu and Kolbuzwe ski's 11144 3 
work on three dimensional ellipsoidal surface deformation 
of a liquid embedded in another liquid, is highly valuable 
for the future work on unstable void growth. Therefore, 
it is clear that the contribution arisirg out of the 
unstable void growth and shear localization between 
cavities leading to ligament ductile fracture should be 
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established. Recent work reveals that the shear localiza- 
tion in a narrow shear band causes concentrated flow 
localization leading to ductile rupture through profuse 
micr ovoid formation and growth within the band * 

The application of continum theory of plasticity to 
such homogeneous and inhomogeneous flow localization 
problems is a major subject of interest (Rudhicki and 
Rice c 145 Asaro and Rice 11146 3 and Yamamoto am ) . 

Experimental work on localized plastic deformation was 
considered bjr Green and Knott [Il25 3? El~Sudani and 
Knott [I1433 j Argon [ri48 3> Backofen [1 149 3) Blackburn 
and Williams Cl50 3j Cottrell and Stokes [IlSl3> Beevers 
and Honeycombe Cl52 3 Price and Kelly 13l53 3» ^ 

theoretical model on flow concentration in two-phase 
materials has been proposed by Cook, Rau and Smith (l9o3* 
According to Orowan, the quadrilateral void formation at 
the centre of a tensile specimen during necking state of 
large plastic deformation is due to high shear localization 
at 45° planes. The necking mode of plastic instability 
C154: and the instability due to abrupt void growth at 
the centre of a tensile specimen are normally considered 
as two independent modes leading to fracture C1533, 

Sih 13 155 3 > Thomason j3 l563> Hancock and Mackenzie [5.5'Q 
and Hill and Hutchinson L. 122 3 considered the details of 
the mechanism involved in these two fracture modes. 



41 


Nevertheless, the basic points put forward by Price and 
Kelly remained that the void formation at the centre of 
the specimen observed by Orowan, was mainly due to high 
shear localization on 45° plane. Pig* 2.5 shows the 
nature of shear localizations in the tensionr-compresjsion 
fie Id for the necked region of a tensile bar 159 3 • 

Figs. 2,6 and 2,7 show some examples of shear localiza- 
tions 

The present work shows, if two prolate spheroidal 
cavity-holes are present in an unis.xial tension field 
inside an incompressible elastic solid, then there will be 
one or at the most two high shear strain localization 
regions. There is a particular angular orientation bet^reen 
the two holes at which large scale plastic yielding on 
the void surface will be visible together with a high 
shear localization in the ligament. This model is then 
extended to a plastically yielded material. Our analysis 
on plastic stress field satisfies all the necessary and 
sufficient conditions for velocity discontinuity and 
plastic instability, put forward by Hill, Kachanov and 
Berg-McGlintock, An approximate analysis based on Henclsy-w 
Prandtl Slip- line field construction reveals that there 
exists a small region bet^vjeen the two holes, in which there 
is a steep shear strain gradient. During an incompressible 
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slow plastic flow, such large shear strain localizations 
and velocity discontinuity in a narrow hand are favoura- 
ble for local plastic instability, and possibly tantamount 
to fracture. 




• GHAPTER III 


STRESS IIirSNSITY FACTOR FOR A CRACK EMANATING FROM A 

SEMI- CIRCULAR EIDGE NOTCH 


The problem of determining the stress distri- 
bution in a semi-infinite plate which contains straight 
radial cracks originating from the boundary of a semi- 
circular cut-out machined notch, and which is subjected 
to a prescribed loading, is one of great practical impor- 
tance. Solution of such problems may provide some theore- 
tical justification for the stress relief used in enginee- 
ring practice. During dynamic service loading, small 
fatigue cracks are normally seen to emanate from the root 
of machined stress concentration regions .Cracks ahead of 
welded joints, key-ways, fillets, under-cuts, roots of 
turbiiB blades and gear teeth are a few such examples. 

Many catastrophic failures in actual service situations 
occur due to the presence of such short cracks emanating 
from the root of radiused machined notches. Mubeen and 
Bandyopadhyay [^52, 176, 177 J experimentally studied the 
fracture behaviour of such geometrical bodies in different 
materials. The most general and complete discussion of 
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the subject would appear to be those due to Bowie 
li'Jig gels worth Newman [l34[3j Nisitani 

Smith, Jerram and Miller Cs^U and Rooke and CartOTight 
[3 96^, Only a few analytical solutions are available 
in the literature so far on this type of problems., lliv; 
present work, firstly, demonstrates the application of 
finite Mellin transform to such type of problems, and 
then obtains an asymptotic analytical expression for 
S.I.F* of crackj after solving a Meinsr-Hopf integral 
equation Dsn. 


3.1 Finite Mellin Transforms of Stresses and Displace- 
ments in a Notched Plate 


In terms of polar coordinates ( p , ([>)? a-s 
shown in Fig. 3.1(a), an Airy stress function x(p, 4>') 
satisfies the well-known biharmonic equation D79n 
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Our unknowi x(p,$) has certain properties on the cut-out. 
rim P = a, which are explained in Appendix A. Following- 
Green and Zerna ClSi^ and Love the stresses 

and displacements may be expressed in terms of the Airy 
stress function X ( P ,<|)) as follows; 


Tangential stress, 


a^x 

ap^ 


45 

(3.2) 

Radial stress, cfpp 

li 

a^x 

1 ^ ^ 

+ p ~ 

(3.3) 

Shear stress, „ . 

’ °p (}) 

_ 1 

ax 

'^2y 

1 a X 

(3.4) 

p2 

a<t) 

p ap a^ 


In the case of generalised plane stress, radial displace- 
ment, Up and tangential displacement, are related to 
Airy stress function as : 


B ( i 
P 


Vl 


+ 


P 


=c 


ap^ 


1 a^x 


-V . 

p^” a(l) 


1 dx 


2 ^^.2 "p 


ap 




...(3.5) 


E 


ap 


1 a^x 

pS 


1 ax 

p ap 


- V 


a^x 

a p' 


...(3.6) 


B( £ — a.^. ^)=2(i+v)[; 




So 


"J 

P 


1 ax 


a<t) 


~,2 Y 

1 a X 

p Jpd^ 


t-irj 

...J 


...(3.7) 


E denotes Yoimg*s modulus and "p is the Poisson*s ratio , 

No-w, if we transform radial distance p to a 
new variable w through a finite Mellin transform, then 
the integrated values of transformed stresses, displace- 
ments and Airy stress function in the finite range are 
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defined from Appendix A, as, 
oo 2w-2 

I X (P, $) p dP , for p 1 a 

a 


oo 

= ! x(p 
o 


2\f^2 


dP for P ^ a = ( 3 * 8 ) 


CO 2w 

I axK(P,<|)) p dP , for P La 


a 


a 2 w 2 ^ 

= X ( P dP + J* P dP for P > a 


2w + 1 


2 w+l 


+ p/b/h (w, <!)) 


( 3 . .^ 9 ) 


OP 2w 

X Odd^Pj^^ ^ J P La 


a 


OP 2 vr 

X Ppp(Pj<l)) P dP , for P > a 


o 


app(v,<l)) (3.10) 


2w 

X cjp^(p,(l)) P dP , for P L S' 


oo 


2w 


X op^ (pj$) P dp , for P > a = (w,<l)) ( 3 . 11 ) 


oo 


Sw 


X Up (p,(I))p dP , for p L 
a / 


2w 


(Up ) a " oo 2 ¥ 

+ X Up (p,(l)) P dP , for P > a 


48 


1 (j) = ^72 
for P a 


r- , sk , , „ / , v2k 

T Cl + 4 2 (-I) ■ {-(k-1) Bg^^Ca/P) 


k— 1 y2j3j4j««* 

'2k+l' 


, . .2k+l 

+ k Cpi^+i (a/p) 


T Cl + B 


R (a/p)^C (3.14) 


11^3 j 4 j 5 j 6 j • ) 


n 


The coefficients are calculateci from Bovxie^s result 
£184 3 who has actually computed the values of Bg^ 
and Cl 2 j 5 .^]_ and tabulated them. Table 3.1 gives the 
numerical values of the coefficients, Now, let us 

transform the boundary conditions and get the following 
relations: 

On |) = 0 plane, the normal and shear stresses are zero 
on the flat edge^ 


X = d X /d^ - 0 , (3.3.5: 

and on 4i = ■ir/2 plane, the shear stress is zero. 

This gives, 

d X /d^ = 0 (3,16) 

The tangential normal stress on plane ^ = tr/2 can be 
written as the sum of the stress from p = a to 1 and 
P = 3. to infinity. Thus, transforming the tangential 
normal stresses bn (j) = 'r/2 plane one gets, 
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2w (2w - 1) X (v, ir/2) = F (w) + o’* (-w) 

, v 2w 


(3 *17 ; 


a 2 tj 1 2-w 

•where, F^ (w) = J* ( 72*^ '^p J" p j'n’/2)P dP 


1 Roa^ 4 a Rsa: 

“ ^ 2w+l 2^-2 2 \,j!-4 


5 

«*!.;» JJ. 


.n 


, , + - - + , . 

2w - (r>-l) 


(3.18) 


oo 


and G“(w) = J* p^^(p , ir/2) dP , p > 1 (3.19) 

1 


On plane ir/2, at p = a, the radial displacement is zeroj 


so ( Ud ) 


= 0 and if -we measure with reference 

P sa ' Y 


to (Ui) } then, the transformed relative displace- 

(p p = a 

ment at the notch root iSj 


+ 1. 21^^-l 

U (w) = J <|)(P , ir/2') P dP , p > a (3.20) 

a 

Thusj following Doran and Buchwald C 27 3> VIeinor-Hopf 
type of integral equation can he constructed satisfying 
all these boundary conditions. Appendix A describes the 
nature of the solution of the differential equation for 
transformed Airy stress function. 

From Appendix A, ^ 

5r(w|^)) = A cos (2w -•1)|) + B sin (2w-l)|) + G cos(2w+l)|) 

+ D sin(2w+l)(l) + X (w,<l)) (3.21) 

c 
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The four unknown constants A) B, C and D are determined 
from Bqns. (3,15), (3.16) and (3,17), and then, suhstitU“ 
ting for transformed displacement to Eqn, (3.20), one 
obtains? 

• 4 * 

» + J. H (w) 

G (w) + P (w) + 2E U (w)... = 0 (3,22) 

H (w) 


where , 


H***(w) 

HT (w) 


w 


sin 2irv/ 


( sin^Tn'r - 4w^) 


Thus, following Dor ah and Buchwald L27 Jj it 
is possible to separate the function H (vr) in two regions 
S'*" and S . The function is analytic in S*’’ when super- 
script (+) sign is written on top of a function, and 
similarly (-) sign is for analytic function on S”, 


Eqn, (3,22) car. also be writte.i in the folio wi.n 

form : 

2B. if (w) H'^Cw) + C'‘'(w) = - G*"(w) ir (w) - C"(w) 

= - IlG^ Xw) + C" (w) 3 = K(vr) 

j'O QO ' 

where, f’*’Cw) h"(w) = g’^ (w) + C" (w) 
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C (vr ) 


= T C 


2w +1 


{H** (w) - i ) } + 


s, 

2^lJ^2 


{H (w: 


R Si 

- H" (1) } + {H" (w) 


♦**' , C* ^ -V 

H , - . -■ -i- 


and 


C***(w) 


Rqa 

O 


^ H (-1/2) 

T C — ™ + 

Sw 1 2*vr*- 2 


H (i: 


. ■ 4 

R/* cl <r, 

4 ; O ^ 

3 2 




H^Cv) can be calculated from Titchmarsh an-5- 

Noble formula t 


T + ioo In HhCz)] 


iT (w) = e 


27ri Y ~ ICC 


dz 


z - w 


for T > Pi3 (v J 


This integral can be evaluated either numerically o:r 
by applying residue theorem in complex plane, cons-''^ 
ring poles and zeroes of H(z) function. 


The stress intensity factor can be obtenned 

from C 27 3 •• 


K-r = vAr vAl-a C ( 0 ) 


:3 .24 j 


R 

-where, G '*'(0) s T Ch“ (- - ) a^ h“ (1) + 

2 2 


+ H* ( - *1 + - h” " f' 


2 




3 
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1 


+ 1 

4H ( i ) 


1.1215 -./^2 


H*” (1), h" (3/2 )j H*"(2) etc. have been calculated from 
the properties of generalised factorial function [^185 3 
In terms of Gaussian products | , 


TT Cl 

m/2 > w 




W u/ is 

) e S 3 

m/2 


/■|a- pflCa- ■:^)a- TTca - ^ ) 






rr w/n _ 

a- 1) B 3 

^n 


(3,25) 


where, are the roots of sin ttw + 2w = 0 and being 
the roots of sin ttw - 2w 0 equations, These roots 
have been calculated by Doran and Buchwald 327 185 J . 

Taking logarithm on both sides, the asymptotic summation 
series converges very slowly, but still the series is 
absolutely convergent (Whittaker and Watson j_ 186 3 )• 
Table 3.1 gives the computed values of IT ( --g“) for 
n - 3,4,5, , , Finally, the stress intensity factor 

for a crack emanating from a semi-circular edge notch 
is given by, , 
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^ _ vn>.'a Cl- S \ (3.26) 

j 4 j 5 j • . 

where, Kq = 1,1215 Tx/Vv/"! and for n = 3,4,5,,.,. 
are calculated and shown in Table 3.1, 

Fig, 3,2 shows the non-dimensional stress 
intensity factor for a crack radially emanating from 
the root of a semi-circular edge notch against the ratio 
of crack size to notch radius. These computed results 
agree favourably with photoelastic experimental observa- 
tions, Hence, the three components of stresses, very 
close to the crack- tip, can be determined from Eqns. 2.1. 



TABLE 3,1 Numerical Values of the Coefficients in 


Eqne , (3 . 14 ) , (3 . 25 ) and (3.26) 


n 


H“ ( 2rJ. ) 

Bn 

-3 

2". 65132^' 

0 ,7314 


4 

-2,60824 

0 .3164 

-0.19571 

5 

2,71264 

0.3149 

0.15193 

6 

-0,35248 

0.2839 

-0.01424 

7 

-0.69060 

0.3932 

-0.03220 

8 

0.27900 

0.1635 

0.00464 

9 

0.28240 

0.2616 

0,00657 

10. 

-0.21728 

0.1544 

-0.00265 

11 

-0.12320 

0.2950 

-0.00259 

12 

0.17000 

0.1062 

0.00117 

13 

0 .04776 

0 . 2262 

0.00064 

14 

0.13464 

0.1021 

0.00075 

15 

-0.00896 

0.2456 

-0,00011 

16 

0.10780 

0.0774 

0.00040 

17 

-0.01216 

0.2025 

-0.00011 

18 

-0.08704 

0.0750 

-0.00027 

19 

0.02376 

0.2162 

0.00020 

20 

0.07128 

0.0604 

0.00016 

21 

-0.02960 

0.1859 

-0.00020 

22 

-0.05880 

0.0586 

-0.00012 

23 

0.03256 

0.1955 

0.00021 

24 

0 . 04840 

0.0486 

■ 0; 00007 

25 

-0.03360 

0.1725 

-0.00017 

26 

-0.04032 

0.0480 

-0.. 00006 

27 

0,03380 

0.1820 

. O.OOOi? 

28 

0.03380 

0.0411 

0.00004 

29 

-0.03304 

0 . 1645 

-0.00014 

30 

-0.02856 

0.0407 

-0.00003 

31 

0.03180 

0.1721 

0.00013 

32 

0.02400 

0.0357 f 

0.00002 

33 

-0.03008 

0.1580 

-O.OOOIl 

34 

-0.02048 

0.0354 

-0.00002 

35 

0.02856 

0.1647 

0.00010 

36 

0.01700 

0.0316 

0,00001 

37 

-0.02736 

0,1533 

-0.00008 

38 

-0,01440 

0.0315 

-0.00001 

39 

0.02584 

0.1596 

0.00008 

40 

0.01216 

0.0288 

0.00001 



CHAPTER IV 


ELASTIC STRESSES AHEAD OF A CRACK-LIK&-NOTCH 

■* It is, of major interest in fracture mechanics to 

correlate the local fracture parameters with the fracture 
toughness of a material. Short cracked beams in metals 
normally fail at a fracture stress comparable to the yield 
strength of the material and at the onset of fracture the 
geometry is very similar to a notch. Analysis on a crack- 
like-notch may give some valuable insight. Combined effects 
of bendirg and tension on crack-tip may produce local strain 
derivative jumps inside the material ahead of a crack-like- 
notch. Thus, it is highly important to know the elastic as 
well as elastic- pi as tic stresses ahead of a crack- like -notch 

4,1 Elastic Stress Field Ahead of a Crack-Like- Notch 

The present section describes a method of determi- 
^ ning the elastic stresses at the root of a crack like notch 

Let us consider the geometry shown in Fig. 4,1, An 
infinitely deep external charpy type notch with a small 
flank angle, < ±q loaded by a constant pressure over an 
unit length on its flat faces. The tip of the notch is 
rounded by radius ^a* . The romded portion having 
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included angle, ir - c.< at the centre, is traction-free. 

Thus, our prohlem is very similar to Ohreimovhs wedging 
problem, but with a finite tip-root radius. This problem 
is also analogous to the bending problem of a 45° charpy- 
notch studied by Griffiths and Omxi [159 3 or a hyperbolic 
external deep notch studied by Neuber Ll9 3 , 

We propose to solve this problem by a simple set of 
imaginaJry cutting and welding operations 

make two straight • , through cuts on the bodj^-j one along 
ODA and the other along OEC. Thus, we are l.eft out w5.th 
two identical 90° wedges ^ symmetrically resting on a 
single TT - oC semi-circular edge-notched wedge plate. In 
order to restore these regions to their original form, 
let us introduce cr and t as surf ace- traction distribu- 
tions applied on the plane (and similarly on EC), cr 

I" 

is the normal stress dist^. ibution and being the sheai’ 

o ' 

stress distribution on these 90 - sectors acting on tha 
faces DA and EC, Equal and opposite forces will act on 
the straight boundaries of the semi-circular notched plate « 
Pig, 4,1 (b) explains the cutting and welding operations. 

Now, if wa rejoin all the three bodies together by welding 
along lines DA and EC, and leave it as original body, then 
some displacement mismatch may be generated along tbese 
straight lines. Now, if we do not allow any such displacement 
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discontinuities to happen on these planes, hy satisfying 

the continuity of displacements and their derivatives 

along welded lines DA and EC, then we are left out with 

the original prohlem, hut \-rith an unlaioivn layer of body 

force artificially produced on the ir - c< notched sector. 

The purpose of formuH-ating this problem in such an 

unconventional manner is to prove later thr?t, when the 

1 ° o 

notch flank angle is in the order of g or 1 or even 
less, there will be a mismatch radial outward pulling 
force generated on crack extension plane giving an account 
of Irwin-Barenblatt crack extension force for Obreimov-Mott 
Cottrell problem. 

For 90°- wedge sector plate, we select D as the 

■ # ■ - 

origin, with polar co-ordinate (r, 0), as shown in Fig .4.1. 
The stresses and displacements of this region can easily 
be calculated by applying a standard Mellin transform 
method in the zero to infinity range, as suggested by 
Tranter j Kbiter and Sternberg [123]]] . 

Let X (r,0) be the Airy stress function for this 
90°-wedge with o and t acting on 0^- plane over the face 
DA and a pressure loading on the face DM, with-T over 0 < r 

‘ O 

< 1, on 90 plan^. Then we have 

/ V '1 T+ jF> ^ . -2p+l 

X (r, 9) - — — j (p,6) r d(2p) 

® 2irl s 

Y - 
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Y+i°® 

-i- ; ■ CaCp) oosCsf-TJe + B(p) sin(^lTG 
2rr± Y-ioo 

+ G.(p)cos C2p+rT0 + D (p ) sin (2p+3.T^ U 
-2p+l 

, r d(2p) (4.1) 


2p~l 

where, r is the usual Mellin transform, multiplying 

factor with 2p-l as transformation parameter on Airy 
stress function, as explained hy Tranter Sneddon 

E 143 and Oherhettinger ClsEI * 


Constants A,B,C,D are obtained from the knowledge of 
stresses on 0 =0® and 90° planes. After solving four, 
equations 2p(2p-l) x (p,o°) = F(p), 2p(2p-l) X (p,7r/2) 

o S 

T ' 

= ^ 2p X *(P}0°) = S(p), and 2p x’ (p,Tr/2) = 0 

2p+l ^ ® 

...(4.2) 

the Constants are 


A = 


2 1 

F(2p+1) (2p-sin irp) + ^ S (2p-l)"sin 2'rp - (2p-l) T sin mp 

4p(2p~l) (4p^- sin^TTp) .. .. 


B 


C = 


D = 


S (gp-l ) (2p+sin^irp) + 2 "F( 2 p+l ) sin 2Trp + 2pT cos irp 

4p(2p--l) (4p^ - sin^irp) 

T sin irp + F(2p + sin^mp) - ^ S siii 2'irp 

4p(4p^~ sin^p) 

2 1 

S (2^+1 ) (2p-sin irp ) - ^ F(2p +1) sin 2irp « 2Ip cos irp 

4p (2p+l ) (4p^ sin^p ) 

, . • (4,3 ) 
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Then, the transformed stresses and displacements, , may he 
written in terms of transformed Airy stress function as, 

00 d Xg 


00 ^ 

^ee = I '^gQ(r,9)r^Pdr = "qq(p, 0) = 2p(2p-l) 

00 ^ d X 

“re = {, ^ ^ = “re^P’®^ = 2p 


and for a generalized, .plane stress situation, 

/I St. 
a X, 


00 


M 2p**l *S ‘ , 1 

Up = / Up(r,0) r dr = «• — H — ^ “ (2p~l) 

o ^ 2pE d9. 

(i+2pv)Xg;] 


00 


2p-i /' :. .. 1 

* 0 ^= ! UgCr,© ) r dr = - 


0 


n 3 «- 

r i.’^s 

an llw iii .a w i w ii l ^ Wi i wui igOtt j “ 

2p(2p+l)E de^ 


+ <1 - 2 p (!-'>) + 4 ( 2 + ^)p^}. 


d Y 

1 (4.4) 

de 

By usir^ standard inversion formula for Mellin transform, 
one can evaluate the stresses and. displacements, provided 
E and S are known* Based pn the Ifellin transform, let us 
keep F(p) and S(p') as ipknoTOi quantities, defined on 
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9=0 i.e, = Tr/2 - c</2 plane, a and t are transfor- 
med as. 


oo 

F(p) = X cr 


oo 


ee 


(r,0°) r^^ dr and S(p) = X ct^q (r,0° ) r^-^dr 


o 


2p, 

(4.S) 


Thus, in terms of two -unknown quantities, F and S, \re 
have teen able to determine all the stresses and 
displacements. 


Let us now take the finite tfellin transform for 

the notched rr - oC wedge, with a radiused cut-out 

machined hole at the vertex. From the knowledge of 

< ■ 

finite hfellin transform mentioned in section 3,1 and 
Appendix A , the ^ Airy stress function is 

given by 


x(p,^) = X + X (P,|)) 

field . ccjre 


(4,6) 


where , 


X = 7 a ( -p -1) 

core ^ a 


2m+l 2m 

m (InP/a) ' d 


• .2 (- 1 ) 
m=o,l,2,3, (2m+l)I 


...(4.7) 


1 Y+PO ■ ■ ^ : 

X (P,^) = — — * X (‘Wjf’). cos(2w -l)^ 

field 2vi 

Y-35)o ^ 

\ • t. ' ^-2w+l“: 

+ C ( w,0) cos (2 w+l)$ 3 P d( 2w) 





61 


The constants are determined from the stresses given 

on ({) = 'rr/2 - c</2 plsine, and vanishing shear stress on 

$ = 0 plane. Therefore,' 

„„ 2W+1 

TBa ■ ■ 

(F ( w) + ”2 ) sin (2 M. +1 ) (v/2 - ^/2 ) + 

+ S (•vl'KSw -1) cos(2'W +1) (•rr/2 - °C/2) 


A( ¥,({)) 


2w (2'W-l) C sin 2w (ir-cc) + 2 w sin (ir*' < ) 


C( w,^) 


2w+l 

T0a 

(FCw) + ) sin (2^ -1) (Tr/2~<=C/2) 

d +1 

+ S(w) cos(2w -l)(v/2-c</2) 


2w C sin 2'w(Tr-®() + 2 ■wsin(Tr-'=( )n 

■ . • (4 • 9 ) 


2 w is the usual Mellin transform parameter for the 
notched wedge portion, with p being measured from the 
centr-? of the hole, and angle {) is measured from the 
symmetric plane. On cutting plane DA, the transformed 
variables p and w are basically ttie same parameter 
provided one separates out the notch-stress 

ir/2 - o(/2) = 0°) =P T from the integral 

transformed function, defined in Appendix A and section 
3.1, Hence, if F and S are knoi'jn, then it is a matter 
of simp Ip algebra to determire the stresses and displa- 
cements for the ir-cC notched wedge portion, using the 
properties of inversion Ifellin transform. 



Equal 5 j.ng the normail aaid tangential displacements 
on plane DA , for both 90° sector and notched wedge, 

one gets, 


Uy (r, 0 °), for 0 ^ V - UpCP j'ff/2 - '=</ 2 ), 

f or a < P < 


and 

U 0 (r, 0 ®)j for “ °</ 2 ), 

for a < P < ^ (4.10) 


Solving these two equations, F and S are determined 
explicitly, and by using inversion formula one can at 
once determine o and t on plane DA or plane EC for 
that matter. If F and S are found by this method, then 
one can obtain the stresses and displacement fields at 
the root of a notch. For a given flank angle and tip^ 
root radius, this general method would give the stresses and 
displacements in an external deep notch. 


Fracture mechanics people only accept the stress 
field near the root of a crack, as independent of Poisson’s 
ratio of the material [320 - 27 3 • During -tSae fracture of 
fi high strength material, the notch opening angle and the 
tip blunting are normally very small. So, we are merely 
interested in incorporating the blunting effect together 
with a small flank angle opening to the present knowledge of 



63 


crack stress field. Therefore, the reqTairement of 
fradtiire mechanics demands for the -understanding of 
the stress near the tip of a crack-like- external deep 
notch, particularly ^-rhen =< is very small and the tip 
root radius ^a* is also s-officiently low* 

If we rule out the pcssihility of any displa- 
cement discontinuity mismatch or j-ump on the plane of 
cutting after reweldirg , then from the continuity of 
stresses, the trsmsformed displacement equations in 
terDis of a common transform variable p = w , would 
give as, 



?)hus, froDi these two equations, F and S can be solved, 
jTor a crack- like-notch geometry* Once F and S are 
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knowrij o and t stresses can te found from the inversion 
theorem of Mellin transform* For the equilihrium of 90° 
sector, all resultant forces and moments acting on the 
body are equated to zero, and thereby p is determined 
numerically for a given value of a and <=< , ^ 

determined from the traction-free boundary condition 
X (a, =0 on the rim, followirg Green^Zerna 

given in Appendix A, Finally, tangential normal stress 



radial normal stress and the shear stress 


are evaluated from the inversion theorem of Mellin 



transform* The evaluation of stresses, involve the 
application of residue theorem, considering the poles 
of the denominator functiors in the complex plane, and 
these are described in Appendix B , 
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4,2 Computed Results : 

The radial, tangential and shear stresses are 
evaluated hy applying the residue theorem in complex plaae. 
provided the roots of the two equations are taaown before 
hand. The first few roots of sin 2z('?r-c<) + 2z sirp( = 0 
equation are given in Table B.2 for notch flank angle = l" , 
Whereas, the roots of {{)(z) - 0 equation for 1° flank an.gle 
are reported in Table B.l. The discussion is confined to 
o< = 1° notch flank angle, tmless otherwise stated. Only the 
misfit stresses are calculated at different notch angles. 

The contours of the integral in all cases were closed over 

■ 'I 

a seml-circiolar path, for which the contribution vanishes. 

The path of integration was chosen in the range of 
1/2 < 2 < 1/2 + X , where X is a fixed value for a 

given notch angle =<, In this strip, the functions remain 
analytic. Hence, by sumndng up all the residues, the stres-ean 
are calculated. Fig. 4.2 shows the radial and tangential 
normal stresses on plane very close to the rim in the 
range of a < P ^ 10 a . Equation (B,8) given in Appendix B 

was used to calculate the stresses in the range of a < P< 2a,. 

'0 , ' ■ 

Airy stress function x is mainly dominant in this region V 

core 

The stresses are finite and behave like an external deep 
hyperbolic notch [^19]]. Pigs. 4.3 to 4.16 show the variar' 
tion of tangential normal stress against P/a on 
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'.different planes in certain ranges. As the stresses are 
varying frcm very small to fairly large values,, this type 
of plot was found to be the most convenient. The stresses 
are plotted at an interval of 10° starting from 0° to 80°. 
Stresses on other planes are also caloi|lated but not presen- 
ted here. Figs. 4.3 to 4.5 show the plot for on 0° 
and 46° planes in three different ranges. Similarly, 

Figs. 4.6 to 4.8 show on and 50° planes, and so on. 

Figs. 4.17 to 4.31 show radial stresses cTpp on 

different planes. Pig. 4.17 and 4.18 show the variation 

of ^ppon 0°, 30° and 60° planes for different ranges of 

p/a. Figs. 4.19, 4,20 and 4.21 describe the mismatch 

radial stress near P/a = 100, on these planes. The radial 

stresses are calculated from Eqn. (B.12) based on residue 

theorem. Figs. 4.22 to 4.25 show the nature of Ppp against 

distance P/a on 20° and 80° planes. Similarly, the 

radial stresses on 45° and 50° planes are also shown in 

Figs* 4,26 to 4,31, All these stresses are drawn for 1° 

notch flank angle. The jump radial stress on 0° plane at 

P = 100 a is calculated with a high degree of accuracy. 

There was no mismatch (or jump) in component of stresses 

anyvdiere, cn the contrary for Ppp componerfe the jump is 

seen to be maximum on 0° plane. So, there is a mismatch 

between the far and near field regions for the radial stress . 

# - 
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The shear stresses' are obtained from Eqn. 

(B,13) again calculating the residues. Figs. 4.32 to 
4.35 show the plot for the shear stress on 20° and 30° 
planes. Figs. 4.36 to 4.39 describe the nature of 
on 45° and 60° planes. Similarly, Op^ on 50° and 70° planes 
are given in Figs. 4.40 to 4.43, and so on. The maximum 
jump in shear stress is seen to occur on 70° plane. The 
tangential stress on the rim P= a is calculated from 
Eq,n, (B.IO) by evaluating the integral. This line integral 
is first reduced to a real integral in the range of 0 to 
+ ^ , and then by applying Simpson’s rule the numerical 
values are obtained at different angles. Fig, 4.47 shows 
the nature of rim stress on p = a for 'the range of 
- 89.5 < ^ < + 89.5 , This agrees favourably with the 

results of others Cl^l - 184 3, Fig. 4.48 shCT\rs the 
different regions of stress fields summarising all the 
results. There are five regions of different stress fields 
ahead of 1° flank angle crack like notch . The stresses 
beyond ^ > 89.5° i.e. on the 90° - plate, are also 
calculated from Eqn. (4,1) for a few limited cases. First 


few typical roots of equation <{)(z) =0 are shown in 
Table B,1 for c( = o.2° and 0.6° notch flank angles. The 
values of the maximira jump stresses are calculated at 0.1, 
0.2 , 0,3 , 0,4 , 0,5°, 0,6°, 0.8°, 0,9 , 1° etc. to study 
the Variation of ndLsmatch, Pig, 4,49 shows these jump 
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radial stres? (on 0 ° plane) as -vrell as the j-ump shear 
stress (on 70 plane ) • The nilsf it energy for the entire 
body is shown in Fig, 4,50. It is interesting to note 
that the decrease in misfit energy for a crack like notch 
has maxima first at 0.2^ and then again at 0.6° notch 
flank angle. 
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DUCTILE FRACTURE DUE TO lOTERAGTION OF 
PROLATE SPHEROIDAL VOIDS 

Normally, in a ductile material, a large plastic 
zone size is developed "below a notch root. By plastic 
constraint there will he some increase in local tensile 
stress ahead of the notch, hut the constraint effect may 
not he so high as to nucleate or to propagate micro- 
cleavage-cracks, rather it could initiate micro-voids 
more easily. The unstable growth of voids and the stabi- 
lity of the ligament region between two adjacent holes 
usually determine the criteria for void coalescence at 
large strains. 

The process of ductile fracture involves, firstly, 
the mechanics of deformation due to surface interactions, 
secondly, shear-localizations in the ligament between two 
holes, and lastly, high local normal- strain-concentration 
on the void surface 55, 132 ][] . The interaction effects 
between two growing spheroidal holes in space are studied 
in relation to the stability of deformation in the ligament 
region. 
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5.1 Shear Localization and Ligament Instabili-ty Between 
Two Prolate Spheroidal Holes in an Incompressible 
Elastic Solid Under Uniaxial Tension 


From the previous discussion it is clear that the 
shear strain localization between two three dimensional 
holes could cause a local instability inside tlie material. 
Therefore, it is of great practical importance to determine 
the maximum shear strain concentration between two sphe~ 
roidal holes. An attempt has been made in this chapter 
to calculate the exterior strain field at different points 
between two spheroidal regions. The method used here is 
fairly general and would be applicable to two ellipsoidal 
regions. Cracks, cylinders and spheroids are merely 
special cases which can be derived from the general 
formulatiof. The numerical computations are confined to 
a pair of prolate spheroidal holes in space inside an 
incompressible elastic solid which is shown in Fig. 1,14. 
McClintock Cl34n _ emphasised the importance of this 
particular geometry in relation to the shear localization 
between two voids leading to ductile rupture . 


I^ an infinite elastic medium, the displacement u^ 
an exterior point for an ellipsoidal region is given 
^ 0-022 : 



1 

16 IT G (1 - v ) 


T 





irrG' 


(5,1) 
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where, <|) = dv / It-IQI is the harmonic Newtonian 

potential, and 'P is the hiharmonic potential of attrac- 
ting matter of unit density filling the ellipsoidal volume. 
Therefore, 'P — l*^**To Idv , the integrals being 

taken over the volume V. p denote stresses in usual 

X J 

tensor notation. A repeated suffix is alx^ays summed over 

the values 1,2,3 co-ordinates and suffixes preceded by a 

coma denote differentiation notation. The superscript 
T T T 

T on u^ > ®i3 > Pij denotes the transformed displa- 

cements, strains and stresses after loading, whereas the 
superscript c on u^ , e^^ p^^. denotes the constrained 
displacements, strains and stresses respectively. G is 
the shear modulus and v being the Poisson ^s ratio of the 
material. The definitions of constrained and transformed 
strains are explained by Eshelby [3 161, 187]^ 


The blharmonlo potential ♦ ogn again be expressed 
in terms of harmonic potentials of solids bounded by the 
ellipsoidal surface having variable densities as. 


^ - 2x^ X 


X,. * dv 


V 


t-t. 


+ X 


r ’^dv 


V lT--n 


(5.2) 




la other words, ^ = (x^ + x ^ - 2 

12 

" - 0 g ) 


Yc 


(6.3) 
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iliy^ = rxgdv/lT-T.l , <>zo = -fy xsdv/ |T-r,l , 

d s = ]• xi^av/ |t-t, 1 , ilva = )■ 1 t-t,1 

Xqy ^oy 

and |) 2 = I 2 : ^ dv / IT -To I 

2© V 

The volume integral is taken within the ellipsoidal region 


X t € • 



with origin 0 


at 


X = 0 , 

1 ’ 


X2 = 0,Xg=0. 

The suhscript notation on |) has the followir^ meaning, e.g, 

$p gives a ha^i^onic potential outside an ellipsoidal region 

having density distribution inside the body varying as p_ 

th 

law, (as commonly used in gravitational field calculations). 
Thus, <{^(or 4) represents the exterior potential at a point 
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having unit density of the ellipsoid, $ stands for 

3£o 

exterior potential with ellipsoid density linearly varying 
with al'^ng axis -1. Similarly, 2 represents 
exterior potential having ellipsoid density increasing 
with x^^ law inside the body along axis - 1, and so on. 

Applying Hookers law on Eqn. (5.4) and then 
simplifying in terms of strain components, the displace- 
ments are expressed as : 


1 r- T T T 

87 r(l-v) °12'^,il2 ®13 ’^,il3 

®21 ’^,i21 ^ ’^,i22 ®23 '*',i23 ‘^ ^31 

^ ^ ^ T . -I ■ i r- T . T" x 

®32 '^,i32 ®33 ^^,133 J " W' L®il'l^,l ®i 2 ^,2 

+ ®i3 ^>,3 H- ^ ®11 ®22 f33 ^ <l^,i 

... (5 . 5 ) 


So far our discussion was confined to potential 
at a point due to a single ellipsoidal region. Ncn-r, let 
us consider two ellipsoidal regions in space, with axes 
1)2,3 of the first ellipsoid being parallel to axes 1*, 2*, 
3* of the second ellipsoid. The potential at a. point P 
in the ligament between two such regions will be the sum 
of the potentials due to the first body and the second 
tody, provided the inside density distributions of both 
the ellipsoids alopg their parallel axes are identical. 
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Hence, at an exterior point P, the displacements can 
be calculated from the sum total of two potentials. If 
all the strains and potentials are calculated with respect 
to origin 0 of the first "body, then the same Eqn.(5.5) 
may be used with as equivalent to 


<t) = (b^ +4) 

1st body second body w.r.t. origin 0 


2 , 2 , 


Similarly, 'll should be replaced by, ^ 


Where, <l)^^ = <|)^ _ . , , 


o 1st body 


o second body w.r.t, 

origin 0 


2 = <1^x2 ■ + 2 

o “^o 1st body o second body w.r.t. origin 0 


and so on. 


The concept of additive potentials of identical 
objects is well known in the gravitational and potential 
field theories [[188 3 • -The application of multipole 
potential expansion theorem in two body problems is 
explained in Appendix C. 

Let us now determine the strains at an exterior 
point P, due to the interaction of two identical prolate 
spheroidal cavities in an Infinite medium with geome- 
trical configm* at ion as shown in Fig. 5,1. The tension 
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\ 


loading field is remotely applied and parsillel to the 

major axes 1 and 1^ Thus, the applied strains are 

■ A A 

along axis - 1, - v along a^cis - 2 and « v along 

axis - 3. For a cavity, Eshelhy»s ellipsoidal inhomoge- 
neity eigen- strain transformation problem simplifies to 
Cl61, 187 3 s 


A 

e^ . — G . . = e . . 

13 10 


(5.6) 


If we confine our calculations only on Xg = 0 plane, 

Eqn. (5.6), after substituting the additive potentials, 
reduces to the following form ; 


1 r- ^ 

®11 “ ^ Z: C 'f'.l 




4v C'f',1111- 3(j)^3_^) + ii) 

"12^ ,1121 \i2il ” ^ ^ ,21^ “ ®33 ^,11-3 

...(5.7) 


1 1 - T - 

Q _ + -. = ^ 


22 


2 4.-JI- 


11 > 


-1 


®12^ ,2122 



-T 

— 

2112 - 4)^22^ 

+ e_ ( 
22^ 

,2222 

«■» 

mm 

- 

, + 'i' «. 

^ , 2212 


-^33^ 


-3<I) ) 


J- 

®12 


L6 Tn i (,1) 


• . • (5 .8 ) 


8V>-"11 3 <1>,1§) + O-.-,-, - kry,}} 




**^ ^82 ^ ^ ’^ ,1222 "" K12) ^ 


■ "if ■ 

1 O ) + 


,2221 


,2111 '*', 21 ' 

si 21) } 


33 ,12 + ®12 ^ f^ ,ll22 + ^.1212 “ 2(|) op ) 
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JL 1 ^ 51' ^ T 

®33 ' 2 ’ ®13 ”^’^23"^ 

T 

where, -eij Poisson’s ratio v was taken as 

1 for an incompressible elastic mediimi. 

2 


Solving Eq^ns • (5,7), (5,8) and (5,9), the three 

unknown strain transformation ratios e^ / 

^ A ^ A 11 11 ’ 

®22 ^ ^11 calculated. From 

Ke Hog’s CissJ result, the expressions for the exterior 

potentials of the first prolate spheroid and that of the 

second prolate spheroid with respect to the origin of the 

first, are given in Appendix C. The potentials with 

different variable densities are calculated for the 


first spheroid and also the second spheroid with respect 
to the origin of the first spheroid. From the knowledge 
of raultipole-potential-expansion lll67, issj , the 
expressions for the asymptotic expansion of the additive 


exterior potentials f), 
^z2 given in 


t > 

o y o 
Appendix C. 
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5.2 Computed Results : 


_ T _ T ■ T 

Firstly, the expressions for e^^ , e^^ and e^^ 

are obtained by solving three Eqns. (5.7), (5.8) and (5.9). 

Secondly, the shear strains e^^ are calculated for '^= t 

different a/b ratios with "known values of ■=( and 6. 

are given in Eq.n. (5.3), whereas etc are menr- 


tioned in Appendix C. Thirdly, tlia numerical differentia- 
tions are carried out for the functions and ¥ in a high 
speed computer. Once these numerical differentiations are 
completectby central difference formula, it would be possible 
to calculate the ligament strains at any point on =0 
plane. These shear strains are computed, and then finally 
plotted in Figs. 5,2 to 5.10. The strains are to be 
measured, from the line PQ by drawing a vertical line 
parallel to the major axes of the holes. Fig. 5.2 shows 
the shear strain build-up on the line PQ between two 
prolate holes with 00* = 1.00, b = 1/40 and prolateness 
ratio a/b = 1.20, for a position of 0 = tan"'^(S/c<) = 35°. 
The vertical scale is drawn non-dimensionally for a unit 
remote applied normal strain. Similarly, the shear strain 
jumps are also calculated for 45°, 60° and 80° for this 
prolateness ratio a/b =,1.20. For a prolateness ratio of 
a/b = 1.60, the shear strains are again calculated for 
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0 = 45°, 50°, 60°, 70° and 80°. The anguL^r position, at 
which the shears strain becomes maxim-um, is also presented 
in Table 5.1. As our ..calculations are confined to a rubber 
like material (v= 1 / 2 ) only, the strain- derivative jumps 
are discussed in the next chapter. 

T 

The normal strain on the hole surface at 

the equatorial plane is also calculated from Equations 

-■t 

(5,7), (5,8) and (5.9). At the same critical angle, the 
void surface normal strain at point P rapidly rises. This 
rise in the void surface (at P ) normal strain on = 0 
plane is shown in Figs. 5.11 to 5.16 for a number of cases. 
As the material is assumed to be a rubber like incompressi- 
ble elastic solid, there will be a sharp notch like 
expansion from the void surface on X 3 = 0, at point P, 
producing an abrupt displacement jump with a very high 
local strain (or stress). The strain on the void surface 
can be calculated from these plots by adding the ordinary 
eigen- stra|.n on the void surface, which is also given in 
Table 5.1. Table 5.1 describes the details of these critical 
aisles 8 Q for various prolateness ratios. At this critical 
configuration, there will be a shear-folding on the ligament 
at point X producing flow*- localization. 
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CHilPTER VI 


DISCUSSION 

The basic conception of fracture involves three 
stages atomic, microscopic and macroscopic. The relative 
importance of each stage obviously depends on the type of 
material and also the nature of the problem. For a truely 
brittle elastic solid, the tip of a very fine crack like 
notch expands to satisfy the equilibrium conditions before 
breaking atomic bonds in a small cohesive zone [193 • 
Whereas, in the case of a high strength elastic-plastic 
small scale yielding material, the notch tip may not expand 
but an avalanche of dislocations movement is seen to occur 
first very close to the elastic-plastic interface inside 
the plastically yielded small zone on o'^-plane, before an 
abrupt notch surface fracture instability 3, 48, 49, 66 
5?:. However, for a ductile material, one gets a large 
global plastic strain by accomodating local discontinuous 
fracture with the nucleation, growth and coalescence of 
voids. The characteristic feature of ductile fracture is 
fibrous, with a series of concave depressions as dimples 
Css, 83, 1493. 
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In the present work, firstly, the advantage of 
•using finite Mellin transform has been demonstrated, which 
is found to be suitable to resolve -the boundary conditions 
more easily for crack and notch problems. An example is 
also given evaluating the stress intensity factor for a 
crack radially emanating from the surface of a semi-circular 
edge notch in a tension plate. Secordly, this method of 
integral transform is applied to satisfy all the boundary 
conditions for an infinitely deep external crack like notch, 
opened under a constant pressure applied on flat faces. 
Elastic stresses, mismatch jump- stresses, misfit energy 
are calculated for small flank angle openings. The review 
of previous literature on an yielding material suggests, 
micro-cracks aid micro- voids formation ahead of such 
geometries before fracture. So, finally, the interaction 
strain field between two prolate spheroidal cavities in 
space has been studied under uniaxial tension field inside 
an incompressible material, to find the conditions for a 
shear- localization in the ligament and an abrupt increase 
in normal strain on the void surface; 

6*1 S.I.F. of a Crack Radially Emanating from semi- 
circular edge notch in a Tension Plate. 

a sharp tip crack emanates from the surfade 
of a machined hole or cut-out, ttie notch stress field due 
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to the geometry of a hole appears to. be essentially govern- 
ing the propagation of the crack 3?]]. This means, 

the stress intensity factor for such a crack is strongly 
dependent on the local notch stress field. The participa- 
tion of notch vanishes when the crack becomes somewhat 
large compared to the notch radius size. Fig. 3.2 shows 
that, when the crack is about seven times the edge notch 
radius, there is no appreciable effect of the cut-out on 
the stress field close to the tip region. The nonr-dimen-* 
sional stress intensity factor always remained less than 
unity, and increases less rapidly compared to Bowie’s 
problem 11323. Figs. 3.3, 3.4, and 3.5 show the 
isochromatic fringe patterns observed near the root of a 
fine razor blade slit- cut emanating from a semi-circular 
edge notch surface. Fig 3.3 confirmsthe predominant 
effect of the notch stress field on isochromatic fringes, 
especially when the slit length is small. This is clearly 
demonstrated showing a change in the shape of isochromatic 
fringes. On the contrary, our stress field pattern is 
almost identical to that of a long crack observed by Wells 
and Post Stresses at the root of the slit was 

measiired with sufficient accuracy by observing high order 
fringes very close to the tip. S.I.F.s calculated theore- 
tically agree favourably with the photoelastic experimental 

measurements. 
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6 2 Elastic Stresses Ahead of a Crack Like Notch 

The determination of elastic as well as elastic- 
plastic stress field near the root of a crack-like-notch 
is very important not only to imder stand the mechanics of 
local deformations, but also to obtain global specimen 
fracture criteria in the presence of crack tip blmting 
and opening. Layers of elastic body forces are being 
generated below the roots of a crack like notch by imaginary 
cutting and welding operations 179ll» No weak dis- 

continuity in the form of strain derivative jump has been 
allowed on the pl.ane of cutting after welding. Furthermore, 
near the unloaded radiused rim, no weak or strong disconti- 
nuity has been permitted, and thereby plane section will 
remain plane near the notch root. 

The present analysis is valid for generalised 

plane stress case, but can certainly be used for plane 

♦ 

strain elastic situation. The influence of outward radial 
mismatch stress on 0°-plane and the mismatch shear stress 
produced on Vo'^plane is examined carefully from the mechanics 
of fracture. 

The combine d actions of bendirg and tension have 
produced compressive stresses and a pure shear field region 
at some distance from the notch tip . On the radiused rim 
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iE sin ('=(/2) component of the load gives a "bending action 
in the presence of tension opening load T cos ('=</2). 

The stresses very close to the rim siorface (i.e. 
in the core region) havV^heen calculated from Airy stress 
f motion ^ ^ an^ a logarithmic nature of stress 

field occurs in the range a < P < 2a. The value of the 

rim stress "^^ogether with its higher derivatives 

must be computed in order to determine the stresses in this 
region. The nature of logarithmic stresses and the stress 
distribution patterns agree favoirrably with the results 
obtained by others C25, 73 3 • There was no difficulty in 
applying the residue theorem in the range beyond p > 6a 
to find the tangential, radial and shear stress components, 
because the first few real roots were found to be sufficient 
enough for an accurate result. The calculations for 1 ° 
flank angle notch with 1% tip root radius bl-unting are 
only presented. Similar calculations based on this 
procedure can easily be carried out on other smaj-l notch 
angles. For higher notch angles, the complex roots of 
equations (^(z) ~ 0 and sin 2z(ir - cc) + 2z sin «=< = 0 
should be considered to find an accurate stress field. The 
accurate determination of the components of stresses in the 
range of a < p < 6a involves the residue contributions 
arising out due to the complex conjugate roots of the above 
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This has been explained in Fig. 4.48 vhich shovs that there 
are five, distinctly different stress field regions ahead gf 
a crack like notch. These are s 

(!) A core region, in -which the stresses 
behave like an external deep hyperbolic notch (Figs. 1.5, 

4,2 and Neubei* Logarithmic nature of stresses 

are seen to occur in this region. 

(ii) Then over a fairly large portion, a crack 
like stress field exists. A rapidly decaying nature of 
stresses is noticed in this region. 

(iii) The radial stress changes its sign and 

becomes compressive between planes 47 i ° ^ ^ < 70° and 

beyond P ^ 100 a, but the, tangential stress still 
continues to be tensile in nature. 

(iv) Beyond ^ > 70° plane and P < 100a, at a 

far field, the signs of these stresses will be reversed. 

This means, in this region, the tangential normal stress - 
becomes compressive, whereas the radial component of the 
stress -will be tensile, 

(v) For ({) > 70° plane and beyond P > 100 a, 

both the tangential and radial normal components of stresses 

ff 

become compressive due to the bendirg action, leaving a 
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pure shear field region over a small portion. 

The nature of ( on the radiused rim 

surface is calculated hy evaluating the integral and 
pre sented in Fig . 4 .47 . 

6,3 Mismatch Stresses and Misfit Energy for a Crack Like 
Notch 

To understand the stability of the geometrical 
configuration for a crack like notch upon loading, it is 
necessary not only to satisfy the equilibrium equations, 
boundaiy conditions, compatibility conditions, finiteness 
of stresses and strains, but also it is highly important to 
stU(^ the finiteness of strain- derivatives on the body 
everywhere, lAJhen the flank openiig angle for a crack like 
notch is very small, it is seen that on the crack extension 
plane all the above conditions are not fulfilled . For a 
brittle elastic solid it is difficult to generate a genera- 
lised Somigliana dislocation on the parent body in the form 
of a strain- derivative ^'ump* So, the mismatch outward 
extra pulling radial stress will expand the notch tip. For 
such materials, the misfit energy released will be required 
to create the new surface as proposed by Obreimov [”6 H i 
Frenkel C 293 and Cottrell , 
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Calculations for mismatch stresses were done at an 
interval of 0.1° notch angle. The roots of the necessary 
equations are calculated first, and then from residue 
theorem, the mismatch stresses are determined for each notch 
angle. Fig* 4.49 (a) shows the variation of radial mismatch 
stress on 0°plane against notch flank angle openings, whereas 
the nature of mismatch shear stress on 70° plane is shown 
in Fig. 4.49 (h). The total misfit energy of the entire hody 
is shown in Fig. 4,50. The mismatch in a crack like notch 

rapidly decreases as the notch angle is increased beyond 8°. 

. o 

The decrease in misfit energy becomes maximum first at 0.2 

o 

and then again at 0,6 ai^le. These are favourable for 
instability. 


If c is the new crack length emanating from the 


notch tip being generated by abruptly releasing the misfit 
energy on 0*^ plane due to mismatch outward radial stress, 


then, the new crack length c is given by 

1 L (c^OD )m 1 

2 Y c = PP M (6,1) 

C JS Ae Cro-CK uvy't't LocucLiyt^ 

where, (Cpp is the misfit radial outward pulling stress, 


E is the Young's modulus and Y being the energy required to 
create a new surf ace per unit areaj/?e>" 44.y\lt 5 jbec/^£>^ 


This additional displacemert for the notch tip on 
0 plane can also be approximately calculated as^ 



89 


U = 100 a - (6.2) 

E 

where a is the notch tip radius . 

As T^a is normally small, Bowie type expression 
for the stress intensity factor can he used to calculate 
the fracture load of a crack like notch, 

K = 1,1215 Ts '.rir -/U (6.3) 

I 

where, s is the geometrical factor, as explained in Chap. Ill, 
considering the stress concentration and the stress distri- 
bution near the rim on 0° plane. 

Hence, it is clear from Eqns. 6.1, 6.2, and 6.3 
that, when the mismatch radial stress on 0° plane reaches a 
critical value, the fracture of a crack like notch could 
occur. For an elastic material this is also equivalent to, 
when the misfit energy reaches a critical value, crack 
extension occurs from the root of a crack like notch. On 
the contrary, for an elastic- plastic material, it is easy to 
produce a strain-derivative jump and generalised Somigliana 
dislocations on the crack extension plane, and the mismatch 
stresses can be responsible for the initiation of fracture 
ahead of a notch. 
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6.4 Plane Strain Elastic-Plastic Stresses and Mechanism 
of Crack Extension at the Root of Yielded tip Crack 
Like Notch 


The elastic stresses for generalised plane stress 
can he extended to determine the plane strain elastic- 
plastic stresses helow the root of a crack like notch. 

The elastic mismatch between the far and near fields 
remains inside the elastic portion and will not be affected 
due to some plastic yielding near the tip root region. Thus, 
the shape and size of the yieldad elastic- plastic boundary 
can be obtained from the elastic stress field. For a 
Von Mises ideally plastic mate^'ial, the yielded boundary 
is shown in Fig. 6,1. These boundaries ahead of a crack 
like notch are obtained from the relation, 


T 

max 




^ Cc 


^4)(}) ~ ^pp ^2 ^ 



(6.4) 


The way this yielded boundary has been formed can only be 
examined from the velocity of plastic flow field. The 
maximum spread is seen to be on 45^ plane for a plane strain 
situation. The problem has been assumed to be statically 
determinate. At the yi^elded elastic- plastic interface all 
the three components of stresses and their derivatives are 
continuous. Two possible cases can arise, either pure 
shear field region is not within the plastic zone but 



91 


inside the elastic portion^ or pure shear region is within 
the plastic zone producing a hinge due to hendiig inside 
the yielded material. The first case is well known Cl3| 

48, SsHj as maximum normal stress elevation occurs on 
0° plane, in this situation, the crack is likely to extend 
on 0° plane only, with small shear stress on 45 plane, 

1/Jhen the misfit elastic radial stress reaches a critical 
value, the elastic misfit energy is abruptly released ahead 
of the thin last vestige elastic wedge. Thus, the misfit 
is expended to generate an avalanche of dislocations, micro- 
cracks or micro- voids on 0° plane inside the plastic zone 
causing liie onset of local fracture. So, there will be a 
strain-derivative jump /9p at the interface of the 

yielded boundary. From equilibrium equations, this is also 
equivalent to a shear-str«s® derivative jump 
across a thin non-deforming layer producing a velocity 
discontinuity inside the plastic zone, on the maximum 
dilatational region. Hence, for a very small scale yielding 
high strength material, when the misfit radial stress on 0° 
plane reaches a critical value, fracture could occur. As 
socm as the misfit is released on 0° plane, the local 
disturbance thus created ahead of the notch Jpelps to cause 
the notch surface instability ClSO, 153, 154 As the 
misfit reaches maximum first at 0.2® and then again at 0.6*^ 
there will be two distinctly different instabilities for a 
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crack like notch, provided tht-^ notch can open to such a 
large , angle. Then, in this case, the first instahility 
at 0.2° will occur locally inside the material, and the 
second one at 0.6° should he on the notch surface. 

For a large scale yielding situation when the 
plastic hinge forms within the plastically d'^^formed region, 
a strain- derivative jump and velocity discontinuity can 
occiir inside the plastic zone. For a plane stress case, 
elastic mismatch shear stress is on 70° plane, nnd a large 
plasticity spread around this plane is commonly noticed. 
Whereas, for a plane strain situation the formation of 
plastic hinge due to bending action can alone give a large 
strain gradient and local velocity discontinuity (Prager 
and Hodge C66 3 ). Fig. 6.2 shows a fully developed slip 
line field pattern when a pure shear region lies within 
plastically yielded zone and compressive stresses in the yie- 
lded boundary# This velocity discontinuity line in parabolic 
- hyperbolic region is taken from Nadai Sokolovsky 

Ces] , Rice and Johnson and Hill C.SoH • This 

boundary agrees favourably w^(,th the experimental observa- 
tions of others n^l, 1213 . Ihe plastically yielded 
boundary shown in Figs. 6.1 and 6.2 also agree with the 
experimental observations of Hahn and Rosenfield 11413, 
hixon and Visser j” l06 ’]3 , Ewing and Richards UlSl^j 
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Gerl^erich £107 "2 y Bandyopadhyay and Miibeen £52^ and 
Singh , Bandyopadhyay and Murthy £53 ]] . The last vestige 
of elastic wedge below the notch root on 0° plane produces 
a parabolic- elliptic field, which is favourable for a local 
plastic instability. In the case of a fully developed 
plastic zone, the presence of compressive stresses on the 
yielded boundary with a small pure shear field inside the 
yielded material could give the possibility of existing 
an elastic core as a last vestige around this pure shear 
region £50, 65 - 67 3 • The characteristics lines, (which 

coincide with the slip lines for a plane strain case), are 
changing first from logarithmic spiral, (hyperbolic field) 
to straight lines (parabolic field), and then finally to 
the boundary of (elliptic field) elastic core* Plastic 
flow can take place from the two sides of the elastic core 
like a flow past over an aerofoil blade. If no elastic core 
is formed in the case of a fully developed plastic zone, 
bendir^ action alone can produce an abrupt strain derivative 
jump inside the plastically yielded medium. 

Normally, the plane strain plasticity stress field 
and the velocity equations are hyperbolic [”5o3 a 

fully developed plastic zme in a large scale yielding 
material, these partial differential equations can have 
bifurcations giving elliptic equations, especially when 
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rt Cl i „ (Sokolovsky r65'l ). Thus, the presence 

^<p<j, PP inax “• 

of compressive stresses on the yielded boundary can cause 
bifurcations and local necking plastic instability inside 
the deformed zone. The second velocity discontinuity line 
on the last vestige elliptic core is obtained purely through 
graphical constructions, by matching the characteristics 
(i.e, slip lines) of three different fields. As the_boundary 
of the elastic core is unknowi, the problem is fairly complex 
and could even be statically indeterminate. Pending such 
an exact analytical solution for stress and velocity fields, 
in the present work through simple graphical constructions 
(magnified to twenty times scale), a, first hand idea of 
this extremely complex problem is given qualitatively. 

Thus, there are two regions where local plastic instability 
could arise : (i) one at 0^ plane below the notch root, 
where a high dilatation occurs, and the other (ii) on 45° 
plane (for a plane strain case), where pure shear with a 
high hydrostatic tension could exist (Figs. 6.3 and 6.4), 
Approximate stress fields on 0° and 45° planes shown in 
Pigs. 6,3 and 6,4 are for a material of yield strength 
= 0,65 T • In the absence of work hardening for a 
Von fases material, the plastic constrain-t.* factor could 
be in the order of 5 or more below a crack like notch. 

These stress plots and slip line constructions are not all 
that accurate, but would certainly give the order and 
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general pattern for an elastic-plastic deformation ahead 
of a crack like' notch in a ductile material. Micro-voids 
can easily he created in a ductile material in the presence 
of a high hydrostatic tension comhinsd with shear stresses 
alorg 45° planes E '^9, 82, 152, 1533 • 

The local plastic constraint factor on 0°plane 
for a crack like notch could he more than 2,7‘lfor an 
ideal situation. Recent experimental observations on sphere I-- 
dised steel £643 show the p.c.f* could he as high as 
3,5. A high local plastic constraint is also observed 
by Singh, Bandyopadhyay and Murthy E^^E while testing 
the fracture strength of short edge cracked beams in steels. 
These experimental observations on 25 mm thick mild steel 
bars have revealed that when the crack length is in the 
order of 3 mm, an abrupt plastic shear flow localization 
begins on 45° plane at some distance from the crack tip, 
whereas for long cracks of more than 6 mm, the abrupt 
local plastic flow was mainly confined to 0° plane. This 
abrupt flow localization at the root of short cracked 
beams observed through photoelastic coating clearly 
demonstrates that a crack extension could arise along 45° 
plane for such highly yielded crack like notch situation 
E 533 . Rice and Sorensen £733 have reported that Begley 
and Landes C723 also made similar obgerva'tions. 
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6.5 Ligament Instability and Shear-Localization between 
two Prolate Spheroidal Voids inside an Incompressible 
Medium 

A simplified model based on the interaction 

of 

strain field between a pair/prolate spheroidal holes is 
presented in the present work. The major axes and the 
equatorial planes of the holes are parallel, and the two 
holes are situated inside an incompressible elastic medium 
with major axis being parallel to the direction of applied 
remote tension loading. Abrupt strain and strain deriva- 
tive jumps are permitted inside this (rubber like) 
incompressible medium, possibly allowing folding in the 
ligament and formation of sharp tip notches emanating 
from the void surfaces. Such theoretical analysis may 
help to understand the mechanism of unstable hole growth 
and void coalescence leading to local ductile mature. 

Figs, 5,2 to 5.10 show the nature of shear- strain 
distributions in the ligament between two such holes inside 
an incompressible (Poisson’s ratio = 1/2 ) material. The 
first cavils on the left hand side is trying to rotate 
clockwise, whereas the second one in anticlockwise direction, 
as if the two major axes can be aligned along line 00’ . 

The total rotational torque in any ( x = constant) plane, 
is always balanced at all time in the ligament. As the 
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major axes of the prolate holes are aligned to the 
direction of applied tension, this is energetically the 
most favourable configuration. So, the holes can never 
rotate. But, in this situation, the principal axes of 
strain can rotate at a point in the ligament to increase 
local shearing resistance 1383- Hence, in our 

case, there has to be a local shear strain build up in the 
ligament. Pigs. 5.2 to 5.10 confirm that, there exists 
a critical angle 9, at -which the rotational effect on the 
void reaches maximum for a given prolateness ratio. Thus, 
when two prolate holes are situated in space in a tension 
loading field, and if the distances between their two 
parallel equatorial and meridian planes are known, then 
there is a particular prolateness ratio at which the shear- 
strain localisation will be maximum, A local abrupt shear 
strain derivative jump in the ligament will be noticed at 
this configuration. As the material has been assumed to 
be incompressible, gt this critical configuration an 
abrupt local normal strain increase will also be visible 
on the Void surface. There will be a tendency of folding 
at a point in the ligament localizing a large shear strain, 
which further causes a normal displacement jump on the void 
surface with a large local normal strain rise (Pigs. 6. 11 
to 5,16). Table 5,1 describes these critical angles for 



98 


different prolateness ratios. Hence, for a given prolate- 
ness ratio, at this critical geometrical configuration, 
the ligament hetween the holes will have a tendency for 
buckling instability. The decrease in elastic strain 
energy is also ma2:imum at this orientation. It is also 
noticed that, there is no shear localization when 0 < 0^^, 
whereas shear can localize in two small regions for 0 >0-5. 

If the above arguments are qualitatively extended 

to a Von Mises ideally plastic material, then it may be 

possible to give a schematic plane strain slip- line field 

on X3 = 0 plane. Fig. 6.5 shows such a slip line field 

between two elliptic holes with a/b ratio =1.6 and 0 = 

0 

60 and for centre to centre distance = 40 b. Fig. 6.6 
describes the nature of plastic constraint and local shear 
stress b\iild up in a small region. The slip lines are 
terminated at X on Xg = 0 plane, shown by dotted lines, 
assuming that there could be a possibility of existing a 
non-deforming layer of material. Although magnified to 
twenty times and drawn graphically, these slip lines and 
stress plots are approximate only. Inside a plastically 
yielded material, such a large abrupt local increase in 
shear strain, with a possibility of a discontinuity in a 
shear strain derivative, may cause velocity d|,scontinuity 
as well as local plastic instability [^5 , 153 3 > which may 
further increase the void surface strain and lateral growth. 



CHAPTER VII 


CONCLUSIONS 

The following conclusions cen be drawn from this 

work ! 

1. Applying Mellin type intagral transform, the stress 
intensity factor for a crack radially emanating from 
the root of a semi-circular edge notch is obtained 
analytically in a tension plate. These theoretical 
solutions agree favourably with the experimental 
observations on model photoelastic specimens. 

2. The complete elastic stress field for an infinitely 

o 

deep external charpy type notch with 1 flank opening 
angle and •] ^ tip root radius blunting is 
presented. It is proved that, there could be five 
distinctly different stress field regions ahead of a 
crack like notch. Compressive stresses and pure shear 
regions are observed due to the combined action of 
bending and tension. Logarithmic nature of stresses is 
seen to occur very close to the notch tip. 

3. When the flank opening angle of a crack like notch is 
very small, an outwaJrd radial mismatch on 0°- plane as 
well as a shearing mismatch on 70® plane will be 
generated ahead of the notch upon loading. 
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4. For an ideal elastic solid, the mismatch on 0° will 
cause the notch tip to expand to sa;tisfy the equl- 
lihri'um conditions on this plane. The decrease in 
misfit energy of a crack- like-notched body has maxima 

n 0 

first at 0.2 and then again at 0.6 notch flank 
opening angles. The misfit energy on 0®-plane will 
he expended to create new surface producing a stable 
crack extension from the notch root. This occurs 
when the radial outward mismatch stress reaches a 
critical value on 0°-plane. 

5. For a small scale yielding elastic-plastic material 
under plane strain situation, the mismatch on 0° plane 
can produce generalised Somigliana dislocation on this 
crack extension plane. This is equivpl'^nt to a shear- 
strain-derivative- jump across a thin undeforraed layer 
inside the plastically yielded zone at high dilata- 
tional region causing a local velocity discontinuity 
and plastic instability. This misfit energy on 0° 
plane will be released abruptly to create an avalanche of 

c dislocations, micro-cracks or micro-voids. Furthermore, 
it is seen that, the maximum normal stress elevation 
below the root of an elastic-plastic crack like notch 
(with = 1° and a = i ^ of the loading distance 

on flat faces), ooiild be more than 2.71, for a nor>* 
strain-hardening Ideally plastic Von- ]y2.ses material. 
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6. Under plane strain deformation in a large scale yield- 
ing ductile material, pure distortional fields could 
te present on 45° planes -within the plastic zone. This 
could produce either an elastic core, or a thin layer 
of imdeformed material inside the plastic zone on 46° 
planes. The presence of compressive stresses on the 
yielded boundary gives bifurcations on plane strain 
plastic i-ty equations, which could further cause local 
velocity discontinuity and plastic instability at some 
region ahead of the notch tip on 45° planes inside the 
plastic zone. Hydrostatic tension combined with shears 
on 45° planes may help to create micro-voids on these 
pure distortional regions. 

7. The interaction strain field in the ligament between 
two prolate spheroidal holes is studied inside an 
incompressible medim. The stability of growing voids 
and shear strain localization in the ligament are 
reported. For a given prolateness ratio of the voids 
there exists a critical relative angular position at 
■vdiich the ligament is unstable. This means, at this 
configuration, a large shear strain will be localized 
at some point in the ligament which helps to build up 
a high local normal strain concentration at some point 
on the cavity surface. 
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Scope for Fut-ure Work : 

It is necessary to carry out further investi- 
o'ations on the elastic-plastic stress and velocity fields 
ahead of a crack like notch. Some future experimental 
work will verify the presence of elastic core inside a 
fully developed plastic zone. Elastic misfit and strain- 
derivative jumps inside a metal should be understood 
through careful experimental observations on different 
materials, To understand the flow- localizations and 
unstable void growth mechanisms insid*^ an actual plastically 
deformed material, more experimental work is needed. Some 
future work will consider these aspects in more details. 

Although more analytical and experimental 
work should have been carried out on plastic deformations, 
to propose criteria for ductile fracture, however, in 
author's opinion, the present theoretical study on a 
crack like notch and spheroidal hole interfactions will 
be useful to understand the mechanisms of fracture in 
metalSo ' 
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FIG.1.6 zones of pure distortion and high 

DILATATION AHEAD OF A BLUNT CRACK- 
TIP [28] 
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-RLAPPING SLIP- LINE FIELDS IN THE LIGAMENT BETWEEN 
PAIR OF CYLINDRICAL HOLES [ 85 ] 
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POSSIBLE REGIONS OF LOCALIZED FLOW AND SUBSEQUENT VOID 
FORMATION AHEAD OF A 45° NOTCHED CRACK- TIP [85] 
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Fig. 2.3 Graok-tip deformation, plasticity spread and 

isochromatic photoelastic coating obsei^ations 
ahead of a short edge- cracked mild steel beam 

CivsH. 
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FIg'2.6 schematic DIAGRAM OF THE FRACTURE PROCESS IN ALUMINUM ALLOY^ 

■ ' SINGLE CRYSTALS AGE HARDENED BY THE FORMATION OF GP ZONES IJSS 






Fig, 3.3 Isochromatic fringe pattern in d-^rk field on 

plane transraisaion photoelastic model specimen 
at the tip of a fine razor blade slit- cut emans 
from the root of a semi-circular edge notch in 
tension plate. ; 

Machined notch radius = 10,50 mm 
Length of the slit cut = 1.10 mm 

Width of the specimen = 75 mm 

Thickness of the specimen = 3.66 mm. 

Tension load = 71.4 Kg. 

Fringe value, f^ =10.76 Kg/cra/fringe 




Fig. 3.4 Isochromatic fringe pattern in light field on 
plane transmission photoelastic model specimei 
at the tip of a fine razor blade slit-cut emai 
ting from the root of a semi-circular edge-noi 
in a tension plate , 

Machined notch radius = 10.50 mm 

Length of the slit-cut = 2,79 mm 

¥idth of the specimen = 75 mm 

Thiclmess of the specimen = 3.66 mm } 

Tension load = 71.4 Kg 

Fringe value, f^ = 10.76 Kg/cm/frir^e 
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Fig. 3 



.5 Isochromatic fringe pattern in dark field, on 
plane transtaission photoelastic model specimon 
at the tip of a fine razor blade slit-cut: emanar- 
ting from the root of a semi— circular edg© no“t:di 
in a tension plate. 


Machined notch radius 
Length of the slit 
Width of iiie specimen 
Thickness of the specimen 
Tension load 
Fringe value f^ 


= 10,50 mm 
= 9.36 mm 

= 75 mm 

=3.66 mm 
= 26.4 

= 10. 76 Kg/ 001/1*37 inge - 




4.2 TANGENTIAL AND RADIAL NORMAL STRESSES NEAR THE NOTCH ROOT REGION 
ON 0° PLANE. FLANK ANGLE, a=r AND TIP ROOT RADIUS ,a = 1 / 100, 





4.3 TANGENTIAL NORMAL STRESS a-4)cj> AHEAD OF A OR ACK- LIKE- NOTCH ON 
AND 45° PLANES. FLANK ANGLE, a =1° AND TIP ROOT RADIUS, a = 1/100. 



0.28 



5 TANGENTIAL NdRMAL STRESS AHEAD OF A CRACK- LIKE “NOTCH ON 

45*^ PLANES, FLANK ANGLE,a=1^.AN^ RADIUS, a =1 /100. 




TiAL NORMAL STRESS 





.9 TANGENTIAL NORMAL STRESS AHEAD OF A CRACK "LIKE “NOTCH ON 
PLANE. FLANK ANGLE , q: = 1° AND TIP ROOT RADIUS ,a=1 / 100 






4.10 TANGENTIAL NORMAL STRESS o^cS AHEAD OF A CRAGK - LIKE - NOTCH ON = 3D 
PLANE, FLANK ANGLE, a = 1° AND TIP ROOT RADlUS,a = 
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G.4.11 TANGENTIAL NORMAL STRESS Oc|>cL. AHEAD OF A CRACK- LIKE -NOTCH ON ci:>^ 
PLANE. FLANK ANGLE, a = r AND TIP ROOT R ADIUS ,a =1 /lOO T ^ 







G:4.12 TANGENTIAL NORMAL STRESS 0 ^ 1 , 4 , AHEAD OF A CRACK- LIKE - NOTCH 
: cp> = 60° AND 70° PLANES. FLANK ANGLE, a=r AND TIP ROOT RADIUS 




FIG. 4,13 TANGENTIAL NORMAL STRESS Ocfxt, AHEAD OF A CRACK - LIKE - NOTCH ON4>=60 
AND 70° PLANES. FLANK ANGLE , a =1° AND TIP ROOT RADIUS a = l/100 
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FIG. 4.15TANGENTIAL NORMAL STRESS AHEAD OF A CRACK-LIKE-NOTCH ON cp = 80 
PLANE. FLANK ANGLE, a =1° AND TIP ROOT RADIUS, a =1/100 
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FIG.4.16 TANGENTIAL NORMAL STRESS cr^ctj AHEAD OF A CRACK- LIKE- NOTCH ON 
PLANE. FLANK ANGLE, CL =1° AND TIP ROOT RADIUS , a =1/100 . 
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FIG. 4,20 RADIAL NORMAL STRESS (Tpp AHEAD OF A CRACK- LIKE-NOTCH ON cp =:30'' PLANE 
FLANK ANGLE, a -1° AND TIP ROOT RADIUS, a = 1 / 100 . 
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FIG.4.24 RADIAL NORMAL STRESS (Tpp AHEAD OF A CRACK-LIKE- NOTCH ON 4> =20° PLANE 
FLANK ANGLE, a =1° AND TIP ROOT RADIUS , a = l/lOO 






FIG 4.25 RADIAL NORMAL STRESS crp.^ AHEAD OF A CRACK- LIKE- NOTCH ON 4> = 80° PLANE 
FLANK ANGLE , a =r AND TIP ROOT RADIUS, a = l /100. ; 










FIG.4.26 RADIAL NORMAL STRESS a>^ AHEAD OF A CRACK- LIKE - NOTCH ON 4> =45° PLANE 
FLANK ANGLE , a = 1° AND TIP ROOT RADIUS , a = 1 / 100^ I 
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FI6 4.28 RADIAL NORMAL STRESS Op^ AHEAD OF A CRACK- LIKE- NOTCH ON 0=4E,"PLANE 
FLANK ANGLE, a =r AND TIP ROOT RADIUS , a = l/l00. 




FIG 4 29 RADIAL NORMAL STRESS 07 .^ AHEAD OF A CRACK - LIKE- NOTCH ON 4 , = SO- PLANE 
FLANK ANGLE ,a = 1° AND TIP ROOT RADIUS, Q = 1 / 100 












'<,>> l *\ ' 





Q o 





163 





iiii 


FIG. 4.35 SHEAR STRESS (rpcj) AHEAD OF A CRACK- LIKE-NOTCH ON 4> = 30° PLANE. 
FLANK ANGLE, a= l°,a= 1/100 





FIG.4.36 SHEAR STRESS cr^<(5 AHEAD OF A CRACK- LIKE- NOTCH ON <f) = 45AND 60 PLANES. 
FLANK ANGLE, a= 1° a =0.01 



FIG.4.37 SHEAR STRESS CTpcf) AHEAD OF A CRACK- LIKE - NOTCH ON 4 4^^ 

PLANES. FLANK a!nGLE. a = 1 ° ,a =1/100 
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IG.4.39 SHEAR STRESS o-pc^ AHEAD OF A CRACK- LIKE- NOTCH ON 0 =60 PLANE. 
FLANK ANGLE , ct = l“,a =0.01 
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FIG 4.42 SHEAR STRESS crpcf) AHEAD OF A CRACK- LIKE - NOTCH ON c}) = 50° PLANE. 
FLANK ANGLE, ct= l°,a = 0.01 
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FIG.5.2 SHEAR STRAIN LOCALIZATION ON LIGAMENT LINE PQ IN X 3 = 0 PLANE BET- 
WEEN TWO PROLATE SPHEROIDAL HOLES UNDER UNIAXIAL TENSION, 
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IG 5,4 SHEAR STRAIN LOCALIZATION ON LIGAMENT LINE PQ IN X3 = 0 PLANE 
TWO PROLATE SPHEROIDAL HOLES UNDER UNIAXIAL TENSION. 









FIG. 5.7 SHEAR STRAIN LOCALIZATION ON LIGAMENT LINE PQ IN X 3 = 0 PL ANE BtTWEEN 
V OVT^'twO prolate SPHEROIDAL HOLES LINDER UNIAXIAL TENSION. ^ 



8 SHEAR STRAIN LOCALIZATION ON LIGAMENT LINE PQ IN X3 - 0 
TWO PROLATE SPHEROIDAL HOLES LINDER UNIAXIAL TENSION 















FIG,5J1':iNTERACT|0N NORMAL. TENSILE;:'STRAIN ,(?|i )[nt AT 
L . ‘ P PER UNJT REMOTELY APPUEJo'-t ENSILE' STRAIN : 

■ plotted against'' void.', POSITION., ANGLE:--' e in degrees 
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crack - LIKE - NOTCH WITH COMPRESSIVE 




FIG 6 3 ELASTIC-PLASTIC STRESS DISTRIBUTION BELOW THE ROOT OF A CRACK-LIKE 
NOTCH ON 0“ PLANE FOR A VON MISES YIELDING IDEALLY PLASTIC MATERIAL 



Stress inside plastic zone 



FIG. 6.4 ELASTIC- PLASTIC STRESS DISTRIBUTION AHEAD OF A CRACK LIKE NOTCH ON 4 
PLANE FOR A VON MISES YIELDING IDEALLY PLASTIC MATERIAL. 
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FIG, 6.6 STRESS ON X3 = 0 PLANE ALONG LINE PQ BETWEEN TWO PROLATE 

SPHEROIDAL HOLES IN AN IDEALLY VON MISES YIELDING PLASTIC MATERIAL 
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APPENDIX A 


application op finite mellin transform in plane problems 

OF elasticity 

Mellin transform in the infinite range, and parti- 
cularly its advantage for vredge region on plane elasticity 
problems are well knoxm. The definition of finite Mellin 
transform and its application to the various fields of 
engineering have also become popular in recent years 
(314, 15, 190 3 • By making a finite Mellin transform, all 
the boundary conditions for a semi-circular edge notch 
cut-out in a semi- infinite plate or in a notched wedge 
can be resolved more easily. We define our finite Mellin 
transform as i 
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Proceeding in the same manner, for n = 2 one gets, 


s+1 

I = - a (df/dr) + (s+1) a®(f) + s(s+l) f 

B, 8 . 


Similar Ijr for I 3 , and so on. 

Thus, the function f is actually transformed in the region 
a < r < j and at r = a, the value of the function plus 
its derivatives are m'itten separately. 

The hiharmonic equation in two-dimensional plane problems 
of elasticity, in r, 0 polar coordinate is, 
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where, f(r,0) is the Airy stress function, r and 0 being 
the usual polar coordinates shown in Fig. 3, 1(a). Let 
’Va'^ be the radius of the semi-circular, cut-out rim with 
0 as its origin on the vertex edge of the plate. Multiply- 
ing Eqn.(A,3) by r dr and integrating in the range of 
a r < 00 j one gets, 
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Making an assumption that the nature of our Airy stress 
function f (r,0) is similar to x Airy stress function 
chosen hy Green and Zerna Cl^^E ^ semi~circular 
edge notch, having property on the traction- free semi- 
circular rim boundary r = a, as, 

(x )„ = (ax /c'r') = 0 , but tangential stress (d^x/ar^) 

a . T— 

exists on the rim. Thus, for a traction-free semi-circular 
type edge notched plate, followii^ Green and Zerna, we 
make an assumption that the nature of the unknown Airy 
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stress function is, 

(f )„ = Qf/ar)^ = (d^f/di^^ =0 (A. 5) 

Hence, by making finite Jfellin transform for traction- 
free semi-circular cut-out edge notched plate, the 
transformed form of the biharmonic equation reduces to 
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+ (s+l) a ( org0)^= 0 (A, 6) 


It is interesting to note th^t the solution of this 
differential equation has a particular integral in terms 
of an unknown rim stress ( pliis a well known 

complementary function CIl63 3 . 

Therefore, the solution of the differential equation gives 
the transformed Airy stress function as 
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Integrating by parts, the particular integral of 
the transformed Airy stress function can be written in an 
asymptotic form and is, 
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where, Oqq = tangential stress on the semi~circular rim 
'r = a, and /dG^. 

8- cl 

This solution will be useful for a wedge with radiused 
cut-out at the vertex. So the included argle could be 
0 < 0 < 2Tr. It is clear that the solution has two distinct 
parts J - the first part is the usual field solution, and 
the second part gives an azimuthal solution for the core 
stress field very close to the rim r = a. If there is no 
cut-out at the vertex of the wedge plate, then the secord 
part does not exist and our result agrees with usual 
solutions., Selecting a suitable strip in the complex 
plane, < y < for which the function is analytic, 
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then from the inversion formula of Mellin transform 
1^14, 16, 383? Airy stress function becomes, 
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Applying Doetsch’s theorem given by Oberhettinger 31^3 
on the second part and separatiig the azimuthal solution, 
one gets the Airy stress function as, 


f(r,0) = (r,9) + ’^^(r,0) 
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where, x^(r,0) = Airy stress function for the field 
region and x = Airy stress function for the core 
region. 


Y+ioo 

x~(r,0) = — - J* , 3 aCs) cos s0 + B(s) sin s0 
2iri Y«ioo 

+ G(s) iejosiS+B)© + D(s)sin(s+2)03r ^ds 


.CA.ll) 



and X (r.O) = Airy stress function for core region 
core 

i.e. very close to the notch rim, 
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Airy stress function x (r.S) is used to determine 

core 

stresses very close to the rim, and x . ,.(^>9) solution 

field 

will he used to find the stresses in the near and far 


field regions. 



APPEHDIX B 


ELASTIC STRESSES AT THE ROOT OF A CRACK- LI KB- NOTCH 


The hasic formulation of the problem is given in 
Chapter IV. After solving Eqns. (4.11) and (4.12), one 
gets the following expressions for 
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-2w-l 

. (P/a) dw 


fe.4) 


The roots of ^)(z) =0, were investigated in the 
complex plane. It is interestiiig to note that for < < 2°, 
the equation has real roots first. In the range of 
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2 < 10 for < 1° ’the real roots of ({)Cz) = 0 

are given in Table B.l. There is always a real root 

_ + X(o() of ^(z) equation, where is a fixed real 

2 

quantity, and X > 0 always. So, the strip of integra- 
tion is chosen in the range of i + £ <y < i +x - so 

2 2 
1 1 

that the line integral is taken along — < "v < — + x , 

2 ^ 

The function is analytic in this strip. Furthermore, it 
is also noticed that along the semi-circular path the 
integration in the complex plane vanishes. Thus, on 
plane (j) = Tr/2 - ®c/2 the normal and shear stresses are, 


= 2 H T p (P/a) 

n=l,2,3 


^ 2 ^ -1 
n 


n 


for a < P < 1 


-2w„ - 1 -2w -1 

= 2 T p X (P/a) ^ + 2 ■ T y„-P ^ 

n=l,2,3 ^ n=l,2,3... 


for p 1. 1 


CB.5) 


and similarly. 


T = 2 

n=l,2,3 


T B Bj^ .(P/a) 


-2wj^ - 1 


for a < P < 1 


= 2 
n=l,2,3 


-2\'j -1 

T p Bn (p/a) +2 


-2w -1 


n=l,2,3 


T 


for p 1 


(B.6) 


^he coefficients , Bn, X^ and Yn 3-^® determined from 
the residue theorem, for n = 1,2,3,4, 


• * » e 
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TABLE B.l First Few Real Roots of <[)(z) = 0,- Re(z) > i 
For Three Different Flank Angles, oc < 1® 


SI. No. 

Roots for flank openirg angles, <=< in degrees 

0,2° 

0,6° 

1.0° 

1 

0,501'^'^ 

0.503 

0.505 

2 

0.502 

0.505 

0.508 

3 

0.999 

0.999 

0.999 

4 

1.001 

1.003 

1.005 

5 

1.503 

1.506 

1.511 

6 

1.503 

1.510 

1.516 

7 

1.999 

1.999 

1.999 

8 

2.002 

2.006 

2.011 

9 

9.503 

9.509 

2.516 

IQ 

2.505 

2 .516 

2.528 

11 

2.999 

2.999 

2.999 

12 

3.003 

3.010 

3.016 

13 

. 3.504 

3.512 

3.521 

14 m 

1 3.507 

3.523 

3.539 

15 

3.999 

3 .999 

3.999 

16 

4.004 

4.013 

4.022 

17 

4.505 

4.515 

4.526 

18 

4.510 

4.530 

4.549 

19 

4.999 

4.999 

4.999 

20 

••5.005 

5.016 

5.027 

21 

5.506 

5.519 

5.531 

22 

5.512 

5.536 

5.561 

23 

5.999 

5.999 

5.999 

24 

6.006 

6.020 

6.033 

25 

6.507 

6.522 

6.537 

26 

6.514 

6.543 

6.573 

27 

6.999 

6 .999 

6.999 

28 

7,007 

7.023 

7.039 

29 

7 .508 

7.525 

7.542 

30 

7.516 

7.550 

7.584 

31 

7.999 

7,999 

7.999 

32 

8.008 

8.026 

8.044 

33 

8.509 

8.528 

8.548 

34 

8.518 

8.557 

8.596 

35 

8.999 

8.999 

8.999 

36 

9.010 

9.030 

9.049 

37 

9.510 

9.532 

9.553 

38 

9.521 

9.563 

9.607 

39 

9.999 

9.999 

9.998 

40 

10.011 

10.033 

10.055 
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The C9nstant 8 is determined from the equilibrim 
condition of the 90°~ wedge as follows t 
Equating the horizontal forces along the direction of 

applied tension Tj 
oo 

T = X. T dP 

SL 


2 n=l,2,3 2 n=i,2,3 % 


9 1 - i 2 

or, 6 = ^ . 

e 


2 E« / Wy^ 

n=l,2,3 


'f:' 

For ■=< = 1° and a = 0,01, 6- comes out to be 595.00 

ap pr oximate ly , 

The functions used in the above expressions are, 


(B.7) 


L(w,°c) = 2 n^os ir w {cos BwCir^) - cos =< } - sin irw 

{2w(3+ 4w^) sin + (12 w^+1) sin 2w(Tr~oC)}3 

K(w,o() = 4w sinTTW £[ (3+ 4w^) sin^irvr - (12\^^+1)23 
- sin 2Trw cos irw . 

M(w,<=() = 2(cos^ 2 w(t- oc) cos^c() + (sin 2w('ir-o() - 2 vt sin >=() 
{ 2w(4w^+ 3) sin °( + (12w^+l) sin 2w(Tr - ©c) X 
N(w,‘=<) = sin 2 tw {cos 2w(v«o() + cos <=< X+ 2w(sin 2w('ir-o() 

~ 2w sin©<) {(3 + 4w^) sin ir w ~ (12 w^+1) } 
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PCwj c<) = 2 sin irw. (cos 2-w (tt-oc) + cos <=() + cos irwCsin 2w 
(tt-oc) - 2 w sin oc) 

2 2 

Q(w,'=<) = sin 2Tr w sin rw + cos to (4w + sin tt w) 

RCW} oC) = sin 2 to (sin 2w('ir-c<) - 2w sin ®c) ~ 2(4w^+sin^TO) 

(cos Pw (tt-oc) + cos =<) 

I (wj °<) = (2w +1) sin (2-w +l)(Tr/2 - '=</2) cos (2w-l)|) 

- (2-w-l) sin (2\’^-l) (ir/2 - c</2) cos(2w+l) 

J (w, !>,=<) = cos (2w+l) (Tr/2 - o</2) cos (2-w^l)(t) 

. - cos (2w -1) (Tr/2 - o(/2) cos (2w+l) 0 

when $) = Tr/2 - oC/2 , J = Oj and I = sin 2w (ir - c() + 2 t,’ 7 sin ■=( 

(l)(o(,w) = C (4w^+ sin^TTw) (sin 2w(Tr-°() + 2w sin ■=() 

2 2 -- 

+ (4w - sin to) (sin 2w (Tr--c<) - 2w sin °c) J . 

C2w(sin 2w(Tr“C<) + 2w sinoC) -f (4w^ +3) - sin^TO 
- (12w’+ 1)} - .(4w^ - sin^TTw) {■2w(4w^ +3) sin °( 

+ (l2w^+l) sin 2w(Tr-<=c ) )]]] - 4 w^ []] sin 2 to 
( sin 2w(Tr-=<) + 2w sin =< ) + 2(4w^'- sin^irw) 

(cos 2W(Tr-c() + cos cC) 3* C2(4tiJ^- sin^TO) 

(cos 2vr(Tr~oc) ~ qos <=<) ~ sin 2 to (sin 2w(Tr-oc) 

+ 2w sin °< ) 3 

when 1 < p/a < 2. x (p j^) Airy stress fimction can 

core 

he used to find the stresses near the core region of a notch. 



For example 


(a 




core 


2 vn 2 


1 


(1 - a^/p2) C (In P/a) ; 
(In p/a)® 


(In P/a)® 





a VI -1 

5“ 

2 


|( 1 ^ a^/P®) C- - 


(inp/a)' 
2! 




li 


(InP/a) ^ j[y (In P/a) vi 

#a 6! J 


" I "Wa" 


(InVa)® J, ii ClnP/a)'^ „ Iv 


21 


^a 


41 


a 


(InVa) 
61 


a VI 

^a " ■**• 


+ - C-Cln P/a) + 

3 . 


(InVa) 


a 

31 '^^a 


(InVa) 


IV 


(In P/a) 


CT Vi 


- ...n 


51 ^'^a 71 "a 

.. . (B.8) 

and similarly, ^pptnd cTp^ can easily "be determined. 
Applying Green and Zerna’s condition, 

X(a, -- Xj.(a,(t)) + X^(a,t)) = 0 (B.9) 


One determines the notch stress on the rim as, 
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C /T 

B, 


P ^Y+ioo ___ 4^2^ sin^Trvr)^ , 

Y - P® (2w+l )(j) (°( 5 ) [jin 2v (7r-=>c ) + 

“ 2w sin o( J" 


? 2 

(4w'- sin irw) N 


(2w+l) |)(o(5w) 

O mm" ^ 

(4¥ -1) C4w -sin TTvx) R 

!(<!)) 






(2w+l) Csin 2w('ir-o() + 2-w sin =<21 


..(B.IO) 


Writing in terras of real quantityj it is possible to 
evaluate this integral h;; Simpson's rule. The tangential, 
radial and shear stresses are determined from the proper- 
ties of Mellin transform inversion as, 

Y 2w 

^ ^ = -1- X ,1— {2(2w+1) (4w^-1) . 

^ vi Y-ioo 2w+l 


2 2 

^ 4w - sin TTvr sin 2-w (ir-'^c) + 2w sin 

w) 

p . O 

_ 4^ - sin TTw sin 2w(Tr-'=<)+2w siip( ^ 

J(^) ^ £ L + k 2 , 


-2^^^-l 

•I ((!))> (P) dw 


2 2 2 

p ^Y+ioo (4vr - sin TTW ) M __ 

vi V.>ioo (2w+l)L§in 2'w(v<-=<)+ 2wsiip0|) (■=<,■«•) 



H 
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. Y+iOO - 

X { 4 w( 2 v+l) r 


( 4 w^- sin^TTw ) P 


Y«. ioo |) (o( j w ) 

(sin 2 v;(Tr“°() + 2 \r sin ©c) Q 


f)(o<,'w) 




(4'W^- sin^Trw)L 


(sin 2 w(Tr-c< )+ 2 -wsin <^)K 


^ |)(cC,w) ( 4 w^-l) ( 4 v^-l) $(=>C,w) 


1 / r, \’"2'W^1 , 

> ( p) dw 


' ***1 . 


p 2 2 2 

^ +i°° ^ (4w--sin that) (4-w -sin ttw) M 

-i°® (4w^-l) ^(=<}w) [3 sin 2w(ir-o() 


TTl 


(4w^- sin^TTw) H 


( 4 w - 1 ) <!)('=(, w) 

( 2 w+l) ( 4 ¥^- sin^TTw) R 


4 - 2 ¥ sin c< H 


3 


<t)(°Cj w) 




+ 


I-, 


(4w^-l ) 3sin 2W (tr-oc) + 2w sin 


-2'w-l 
} (P/a) dY 7 


... (B. 13 ) 

In stands for d^I /d^)^ and for d^ J/d(t)^ in the 

above expressions. The roots of the equation sin 2 w(Ti>^).. 
+ 2 ¥ sin c< ?? 0 are also determined aid given in Table B, 2 . 
The line integral was chosen in such a manner that , 

. ^ V 3 . ^ \fcr^ 1 where i + x(c() is, the first real 
t of the equation |)(z) = 0 for z > 1/2 . Therefore} 
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at any given angle ([) , ty applying residue theorem the 
stresses can he calculated from the above expressions. 

TABLE B.2 First Few Real Roots of sin 2w (ir-'=<) + 

2 w sin cc = 0 5 Re (xO > 1/2 for Notch 

opening Angle, = 1 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


Roots 


0.5056 

0.9999 

1.5169 

1.9999 
2.5281 

2.9999 
3.5394 

3.9999 
4.5507 
4.9998 
5.5620, 
5.9997 
6.5734 
6.9996 
7.5848 
7 .9994 
8.5963 
8.9991 
9.6078 
9.9988 

10.6194 

10 .9984 

11.6312 

11.9979 

12.6430 

12.9973 

13.6549 

13 .9965 

14 .6670 

14.9956 

15 .6793 

15 .9946 

16 .6917 

16.9934 

17.7044 



APPENDIX G 


MULTIPOLE STRUCTUEE ASYMPTOTIC EXPANSION FOR EXTERIOR 
POTENTIAL FIELD OF AN ELLIPSOIDAL REGION 

’'Multipole expansion" method is well established 
in gravitational, electromagnetic and nuclear charge field 
theories. Recently, considerable attention has been 
focussed on finding out the Boussinesq'- Neuber- Papkovitch 
displacemert potentials of spherical and, ellipsoidal 
regions at an external point by applying ’multipole struc- 
ture’ field theory. Miyamoto dlCS^, Hashin 3^'^^ 

Chen H made successful attempts to solve spherical 

region problems by using multipole asymptotic expansion for 
a row of spheres. For a three-dimensional closed surface 
occupying a region inside an infinite continuum can be 
regarded as the sum of a large number of components of 
equivalent potentials at the origin. Thus, the closed body 
like ellipsoid can be replaced by a series of strengths-, 
monopole, dipole, quadrupole, octupole etc, charges at the 
origin. As Newtonian potentials are always additive at an 
exterior point, the potential due to identical density 
distribution bodies with these multipole-components can be 
added at that point for both the first and second bodies. 



^ » ,rith respect to the origin of the first body 0, 

Therefore, with reopecT, ou 

a ^ultlpole e.peneion is ohtained for each potentxa ^ 

,,, 5 1 shows an iUnstration of »nltipole charge 

tntion for an ellipsoid ooncarnxng the Taylo 

in thr e o <3. imen s i ons 
1 


Let - 
R 


o 5y 


= z (-if (i/nO fc'o ^ 

ir-Tol 

. (l/r) 


£ 
o dz 


+ z, 


n 


(C.l) 


■ nvtprior potential for an 
^ n/v) is the Newtonian exteri P 

where , J „ t < 3^1 ace equation 

pined hy solving Laplace 
ellipsoidal region, ohtainea y 

in three dimensions. 

Kellog Ilf or an 

field potentials for 

~ 

ellipsoid as je ^ ,^^eroid of a > b( = c), 

“—'.r' “ 




t* or a. - 

tr Lc per unit vol^e of the cavity 

Cl/r) . -3- l^^g.3o|-4(aW) 
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where, s = + /(x^ + x^^ ) + x^ 2 + + a^- h^ 

+ 2x^ /(a^- h^) -} + -Ax 2 + Xg + x^ + h^ 

- 2x^ 'f (a^- ) } 

For the second cavity with respect to the origin of the 
first cavity, 0 , the potential (l/r)= , x^, x„) 

per unit volume can be written as 

^ = |)(x^ - cc, Xg - p, Xg ) (C.3) 

Thus, both the origins of the spheroids are situated on 
X 3 = 0 plane, but their equatorial and meridian planes are 
at oc and 8 distances apart respectively (Fig. 1,14). 


As we are given some distribution of charge inside 
a finite spheroidal region, and that we have represent ed_ 
by a density function p (x^^, y^^, z^) of position inside the 
region, then the potential at an exterior point due to this 
charge distribution can be expressed as 


0 _ = 


in (l/B) P( x^, y^, Zq) dx^ dy^ dz^ 


...(C.4) 

Applying multinomial theorem on Eqn, (G.l), we get. 


1 

R 


vn J k 


= 222 


(-1)'"" xj yo z 


o 


,n 


n 3 k (n- 3 *»k)I 3 1 k! 


5x3* dy^ dz^"^"^ 


■ (1/r ) 


. . ,(0.5) 
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Xhusj ¥8 have "been able bo replace tbe total 
potential by a series of niiixtipole asymptotic equival®^"^ 
moment structures. ' The details of the method of evaluating 
dipole, quadrupole, octupole moment tensor triple integrals 
over spheroidal region 8j?e given by Band 
Landau and Lifshitz Cl9l[3* Thus, in our problem, this 
method has been used to obtain' additive potential expre- 
ssions for 3 K 1 K > 5 ^2 end } 2 • 

■*•0 ^O ^0 -^O ^o 

moment coefficients are determined with respect to the 
origin of the first cavity. Hence, the final expressions 
obtained are as follows: 


Vv = <t) + |(a^/5 . <1)^11 2b^/^ . (j) gg) + J: [13/35. 

+ 12/35 .b'^ij)^ 2222 + 12/35. - 

(oc + I &^fJ,ll+ %12^ 

“ I C ( 3/5, oca^) + C3/5.3b^ + S )\222 

+ 3/5.c(b^( n^ 3 _ 2 g + ^, 133 ) + 3/5.e(a^fi "^^^,233^3 


+ 1/24. n 3/35. a^ n^3^3_3_3_ + 3/35. b'^( ^,2222 ^^,2233 

■^^,3333^ ^ ^ ^^’1122 ■^^,2233^ 



4 

+ P n 



+ s + 6 


ni33 


+ 


,2233 + 4 C I a^b(p + o('^8) q 


3 ^2 


+ 4( I b o<p + ocp*^) a 


+ iioTR >^^o 1 r / 3_ 4 


6 ^ “ ,1333 + ^ ocp b « ^,233;] ^ ^ 


+ 2=( a + c{^) ( I e^4 ^ g^2^3 ^ 


22222 


+ 5 C cC|3" + 6 jj 


~-^u . .cp p^.,S ^^,i2222-.5(|-a% ^ 


6 .2^ 2 . . 4 


3 4 


+ - «<-6a- + pc<-) . + I eb’ a _23333 + ^ 


+ 10(|-aV8 + ril + 2c<V+»cV)a 


11222 


+ 10 ( — + H 




+ 10 ( |g b^p + — )<! _2 p233 + 10( 2- ^~ 

+ 2 o(p V + =(■=•6^) n ^^33 + 30( n 


3^2 

< b 


.11233 


-6 

4- O Cd ) 


...( 0 . 6 ) 
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a 

5 


■ 6 ^ 35 


i C c I ■=<a 2 + c< 3 ) n + (^p 2 ^ ^| 2 j ^ 

, 12 l 


<= cTd 


5 ^ > 33 ^ ^ ■*■ P a^/5 ) n 


C a^ n ,-,T + si%^i 


6 ^ 35 ^ 


35 


,122 ^,133) 


" ( I =<®a 2 + « 4 ) 0 + ( I .(p ^2 ^ ^p 3 , „ ^ 


+ 3 (c<3s + I- a2) ( 42 /g ^ p^S 8 2, 


5 - -^118 - ov /a -r ^ 


4* (o( B 9^ Q 


,3,3 + i- Zcic < a ^ 


5 ,233 • 5 ' " -,133-. ■ ^ 7 

+ + /)!J + ( |g. C( + I c< pS® + =< p«)s 

4 

4 3 d ( Id fi j 3333 

^ 0/ ^ — 


>2222 


+ 6( ^ c( aS b2>+o(^^/5 + - o( + c(^p^) Q 


35 

+ 6 ( ~ o( + 

,35 




5 “ ,1133 


+ 4 ( 1 ^ paS^ + I c< 28 b 2 


+ a^p^/5 + cc^p^) n ^^333 + 4( 1^ §a^ + | o<^Ba.2+ c(%)q 

+ 6( c(bV35 + o(§%^/5 j2233 a^VsS 


,1112 


+ oc^P b^/ 5 )fi ^233 


n + OCd" ) 


CC.7) 
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6 ^ is ^ ^’222 + SH ^112^ 

^2 

Bfi- C ( i + 8 ^) n^gll + i dC a ^/5 

(p^+ I pb^) n + 6b^/5. Q ^33 + 3=( ( B^+ b^/5)fJ^y23 


={^)B f 2 -, ^ 


+ (p^+ I pb^) n + 6b^/5. Q ^33 + 3=( ( p^+ b^/5)fJ^y23 

- i EP( ^^ + !■=< + C gg b^+ I sV + p‘^)fl,22g 

+ 3=( (p^ + I Pb^) n ^^22 + 3( a^VaS + cc^'b^/5 

+ p\^/5 + oC%^) a ^^^2 + 3( bVsS + P^^/6)^^^233 

^ . 9 1 4 6 ■ 2 P 4- 


3 2 -1 If-, 3 4 6 2 2 4 

+ i ocpb'^fi -,oq3 + — C ( ^ + - c< p a' + oc 

5 24 35 5 

+ ( I + 2 p\®+ p 5 ) „ ^ |_ 3 

+ 6 ( a^ 6^/5 + — a.%% + | o(^pb^ + c<^p^) q 

35 • o 

+ 68 ( oc ^^/5 + ^ a%^) SJ + 6 ( eb^+ 8 ^^/' 

35 jllSp 35 

^ } 2233 ~ a~l) + “ o(P^a^ + °( b^/5 + o( j 

22333 + p)n ^^222 


2 2 -^ 
: p a + « 


q 2 2 ‘3420 

~ a b ) SJ + 6( -S Bb + B'^bVS) 

35 jllSp 35 

L c< + I =<p®a® . cc^t^/S + =<®eb 


^ + ocS^S + o(" 


22333 ^ 35 ^ 5 

, / 12 . , 12 - Q rv 

35 6 ‘ A 23 o 


^ + 8 ^) n 

5 ,1333 

^,1233 3 + 0 (d ^) 


(G.8) 



C-8 




7 


5 ^,3 


C <1),333 + 


jll3' 


- 5^3- : 


•H* ( 


+ ~ 

24 


1 r- 4 

“ C — “b n 

6 35 


S 2 
5 ‘Sl23, 


,333 ^ 


12 

35 


' M .U 

S” 35 ^ ^ , 

223 

+ 

, l^b 
^ 5 

12 4 
C ( S b + 

4 

„3 ' 
P b 

§ 

) Q 

3 

35 

5 


^ ^ ^,2333 

+ 

i£ 

35 

4 

cc b 


2 2 2 2 


35 " ^ ^, 1133 - 


+ ~ ) Q ^,^^2 


2223 


,1333 


+ 4( 3_=c 2f£^) 


35 


,1113 


, f 12 4 

+ ( c< -b + 

35 


+ I2p ( _ + 

35 


12 2 2 
^ « P b ) R ,1223 


Sb" ocV 


^ P 2 2 

^ ^ + o( p ) n ^ 1123 ^ 


•*0 

+ o(d ) 


0 


(C.9) 
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ix ) * i [: fr . 2 

6 2 35 >11 Qi 


I a%3 (j)^222] 


+ (0(2+ ^ o(a2 + 0(3) Q ^ 

5 5 ,1 

+ CPa/5+o(p)jj T + ir"( 3 4 

>2-' 2 *- 35 

+ I + =<^) (1 + (=<202 + a2b2/35 + ct^b^/S 

+ p^a /5)fi 22 + ( a^V35 + o(,%^/5 ) sj 

’■ ' ,>33 

+ (|o(p a2+ 2o(3e)o^^g3 . i£( 3^4^ 

+ 2a + c<^) (2 + ( ^2^2^ ^ a^P^/5 


3 ^2q ,2 _ 2 o 


+ f b“ ,2 sB + 3< si + I =<2pa‘ 

o 

+ o( p) (5 + 3( ^ o( a^^+ c(V/5 + I o(p^a^ 

3 ^ 

+ ■=< P )^ -loo + 3(@a^ b2/3S + __. 


,122. 3(@a^ b2/35 + ^ ^233 


2., 2 


+ 3( Igct aV +=<V/5)0 „,3+ 0(d"®) 


.. (G.IO’ 


ir 2 ^2. 


To'^ b 


I ^ i 35 ^11 36 ■^'^,223 + ( I + P^) ^ 

- £^(. b^/5 + B^) fi + ( 3 ^^^2 ^ -r 

5 -L g ^ r-* 

+ i C ( aV/35 + =(V/5 + 8^/5 -,. 

+ ( L- b^ + £ + P^) R + ( Sl + 

35 5 V >22 ^ 35 5 ^ ,33 


4 d 
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f I P ,12^ - iCcSccpS ^2^^338 

6 5 

. |_^ . A^/S) . ( 3 3,4 , ,33^2, 36,^ 


^SS2 H- 3( |_ a^b^S + I «Ssfc8 ^ ^ ^^^ 3 ^ 


,112 


+ 3 =( 


( 2 - •b'^ + 6 2 g 4 


36 


5 ^ " ^^122 


+ 3C |g p + Q%^/s) jj 


233 


+ 3 =< ( bV 35 + p^^/ 5 ). 


,133 3 


(C.ll) 


12 

^o 

”"V 


bl 1 r- 

= - «{) + ± r d) 4. 4 ^ 4 . _ v 2 

5 2 35 >.ll ^ ^l^jSsU + — 

5 


“ fn.gH , iC ( a ^^ 3 S 4 ^ V / s ) 


Q 


?11 


+ ( bV 3 J + p ^ b ^/ 5 )a _ + i . b ^ 


> 22 " l 33 ^ 


£-cp b 2 n _ 

5 >12 

3 „,4 . _3.2 


3 - iC ( | o < a 3 b ®. AS / 5 ) n ^^^^ 


+ ( 


Q 


,112 


35 + P "b /6)fi + 3B( o(2b2/5 + sS^/ss) 

+ 3=( CbV35 + ^%^/5 )q 


,122 fiP ‘^^^,233 
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